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ABSTRACT 

Aims. With the aim of identifying rest-frame UV and optical properties, typical of Lycr emitters (LAEs, galaxies with EW(Lycr) > 20 
A) at both low and high redshift, we investigated the morphological properties of the LARS galaxies, in particular the ones that exhibit 
intense Lya radiation. 

Methods. We measured sizes and morphological parameters in the continuum, Lycr, and Her images. We studied morphology by using 
the Gini coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated LARS galaxies at z ~ 2 and 5.7, performing 
the same morphological measurements. We also investigated the detectability of LARS galaxies in current deep field observations. 
The subsample of LAEs within LARS (LARS-LAEs) was stacked to provide a comparison to stacking studies performed at high 
redshift. 

Results. LARS galaxies have continuum size, stellar mass, and rest-frame absolute magnitude typical of Lyman break analogues in 
the local Universe and also similar to 2 < z < 3 star-forming galaxies and massive LAEs. LARS optical morphology is consistent 
with the one of merging systems, and irregular or starburst galaxies. For the first time we quantify the morphology in Lya images: 
even if a variety of intrinsic conditions of the interstellar medium can favour the escape of Lycr photons, LARS-LAEs appear small 
in the continuum, and their Lya is compact. LARS galaxies tend to be more extended in Lya than in the rest-frame UV. It means that 
Lycr photons escape by forming haloes around HII regions of LARS galaxies. 

Conclusions. The stack of LARS-LAE Lya images is peaked in the centre, indicating that the conditions, which make a galaxy an 
LAE, tend to produce a concentrated surface brightness profile. On the other hand, the stack of all LARS galaxies is shallower and 
more extended. This can be caused by the variety of dust and HI amount and distribution, which produces a more complex, patchy, 
and extended profile, like the one observed for Lyman break galaxies that can contribute to the stack. We cannot identify a single 
morphological property that controls whether a galaxy emits a net positive Lya flux. However, the LARS-LAEs have continuum 
properties consistent with merging systems. 

Key words, techniques: imaging - galaxies: star formation - galaxies: starburst 


1. Introduction 

Originating mainly in recombining gas being ionized by hot O 
and B stars, Lyman alpha (Lya) radiation has prove d an excellent 
probe of star-forming galaxies at both low (e.g., [Cowie et al.| 
|2Q11[ ) and high (e.g., |Ouchi et al.|2()T0| ) redshift. 

Send offprint requests to: Lucia Guaita, e-mail: lguai@astro .su.se 
^ Based on observations made with the NASA/ESA Hubble Space 
Telescope. These observations are associated with programme 12310. 


Hayes et al. ( 2005 20091 and jOstlin et al. ( 2009| ) devel¬ 
oped a method to separate the rest-frame UV and Lya emission 
in Hubble Space Telescope (HST) data. In these papers it was 
demonstrated observationally that, in nearby galaxies (z < 0.1), 
Lya emission extends away from the star-forming regions where 
the Lya photons were originally generated, forming the so- 
called Lya haloes. 

Y oung starburst galaxies are expected to be very bright in 
Lya ([Rartridge & Reeble^|1967|). For the past 15 years, star- 
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forming galaxies have been successfully detected at z > 2 by 
identifying their strong Lya emission line. The principal method 


used is the narrow-band technique (e.g., [Cowie & Hu 

1998 

Rhoads et al.[[2000t [Quchi et al.[[2008t [Gronwall et al. 

2007 

Nilsson et al.|2009[): Lya emitters (LAEs) present an excess in a 


narrow band (covering the redshifted Lya wavelength) with re¬ 
spect to a broad-band filter (covering the rest-frame UV contin¬ 
uum). Because Lya photons are easily absorbed by dust grains 
and are scattered by neutral hydrogen (HI), LAEs were thought 
to be a special population of galaxies with special dust and HI 
amounts and distribution. Although extensive studies have been 
carried out to characterize LAE physical properties and their 


special conditions ([Nilsson et al.[[2011 Acquaviva et al. 

2012 

[McLinden et al. 2014[ [Vargas et al.[[2014t [Hagen et al. 

2014 

among the most recent ones), the results have been inconc 

iusive. 


The mechanisms (e.g., interstellar medium geometry and kine¬ 
matics) controlling the escape of Lya photon are still debated. 


The morphology of the rest-frame UV and optical continua 
provides information about galaxy formation an d evolution (e.g., 
star-forming region distribution, merger events, Conselice|2QQ3[ 
|Lotz et aL]|2004| ). Erom the ground LAEs were observed to be 
compact in the rest-frame UV, but multiple components were 
identified in deep HST-resolution images ( [Bond et~ar]|2009[ 
|2012| ). There have also been a few attempts to quantify the mor¬ 
phology of the Lyof emission itself. [Bond et ^ ( |2010| ) explored 
a sample of seven LAEs placed at z - 3.1 (observed-frame 
T(LyQf) ~ 5000 A) by using HST WEPC2 (Wide Eield Camera2) 
E502N narrow-band imaging. They found that, for one source, 
Lyof emission extended till 1.5 kpc (< 1 kpc for the other six), 
a just-slightly-larger scale than the UV continuum. Also, this 
source was composed of two main clumps both in the rest-frame 
UV and in Ly(r. |Einkelstein et 2011a] ) spatially resolved three 
spectroscopically confirmed LAEs placed at z - 4.5 (observed- 
frame A(Lya) ~ 6570 A) by using the HST ACS (Advanced 
Camera for Surveys) E658N narrow band. Two out of the three 
systems showed Lya emission significantly more extended than 
the UV continuum. 


Recently, evidence of extended Lya emission was found in 
the stack of a large sample of Lyman break galaxies, (|Steidel| 
eFaLlIMn] ), generally more massive and dustier than LAEs, 


and of galaxies locate d in overdense and not-overdense regions 
( [Matsuda et al.|2012 ). By stacking hundreds of z - 2.2, z - 3.1, 
z - 3.7, and z - 5.7 LAEs from deep ground-based imaging, 
Momose et al.| ( 20T^ discovered extended Lya emission, with 


scale lengths of 5 -10 kpc. Ho wever, by stacking their sample of 
LAEs [FeIdmeier et~aL] ( |2013| ) just found a marginal detection at 
z ~ 3.1 and a non-detection (z - 2.07). It is clear that depth and 
image resolution were the main factors affecting their results. 

Instead, giant Lya nebulae, powered by active galactic nu¬ 
clei, have been studied by a few authors to assess the role of 


HI scattering and Lya radiative t ransfer effects (e.g., [Humphrey 
jet al.|2013a[ [Prescott et al.|2014| ) 

Local starbursts (p verzier e t al.||2008| [ 2009[ [2 pi0[ [Hayes[ 


et al. 


2013 


2014 


Petty et al.|2014[ Qstlin et al.|2014| ) are unique 


laboratories to study the rest-frame UV in detail and optical light 
distribution, morphology, and to investig ate the mechanisms, 
which allow Lya photons to escape. In Qstlin et al. ( 2014 ) 
(hereafter Paper I) we presented the Lyman alpha reference 
sample (LARS), which is composed of 14 star-forming galaxies 
at z < 0.2. These galaxies were observed during HST cycle 18 
(PI. G. Qstlin) in a set of rest-frame UV (ACS/SBC P125LP, 
P140LP, P150LP) and optical (WPC3/UVIS P336W/P390W, 
P438W/P475W, P775W/P850LP, P502N, P656N, and 


ACS/WEC P502N/P505N/P55 IN, P656N/P716N/P782N) 
filters. Lya maps were generated by estimating the continuum 
at rest-frame A(Lya) = 1216 A, through modelling the galaxy 
spectrum as a composite population of young stars, old stars, 
and ne bular gas. LARS images were published in 
(2013[) (hereafter Paper 0) and further analysed in 


Hayes et al. 


Hayes et al. 


(2014) (hereafter Paper II). We found that the Lya emission pro¬ 
file appeared different from the rest-frame UV and it flattened 
on scales larger than the rest-frame UV. The majority of the 14 
galaxies showed a negative Lya equivalent width at small radii 
and then an increase farther out. We concluded that this was 
due to scattering on neutral hydrogen, which is able to shape 
the Lya emission into the form of haloes. Also, by comparing 
LARS Lya with global physical properties, it appeared that the 
Lya photon escape was favoured in the system with weaker dust 
reddening and small stellar mass. 

The neutral hydrogen content of LARS galaxies was pre¬ 
sented in [Pardy et~ar] ( [2014[ ) (hereafter Paper III). The spec¬ 
tral lines of HI were detected in 11 of the 14 observed LARS 
galaxies and it was also found that the Lya escape was favoured 
in low Hl-mass systems. LARS interstellar medium kinematics 
will be presented in Rivera-Thorsen & Hayes ( [2015| ) (submit¬ 
ted), in Duval et al. (in prep), and in Qrlitova et al. (in prep). 

In this paper, number IV of the series, we address the ques¬ 
tion whether specific galaxy morphological properties could be 
related to the escape of Lya photons and escape in haloes. Note 
that in Paper I the present paper was termed paper 7, due to a pre¬ 
vious numbering. In Sec.[^ we briefly explain how we measured 
morphological parameters (details are give in Appendix A) and 
the process adopted to simulate how some local (z < 0.2) galax¬ 
ies would appear at high redshift (z > 2). In Sec.[^ we describe 
the morphological properties of the sample of local galaxies and 
compare them with local-Universe and high-z galaxy popula¬ 
tions. In Sec. 1^ we study the morphological properties of the 
high-z-simulated galaxies. In Sec. we present the stacking of 
the high-z-simulated sample and compare with high-z stacks in 
the literature. In Sec. and [^ we discuss and summarize the 
main results of the paper. 

Throughout we adopt AB magnitudes and assume a ACDM 
cosmology of {Hq, Qm, ^a) = (70 km s“^ Mpc“\ 0.3, 0.7) as in 
Hayes etal. 2013,2014. 


2. Method 

We present the morphology of the local (0.03 < z < 0.2) LARS 
galaxies and investigate how it would change if the same galax¬ 
ies were observed at high redshift. As explained above, in Paper 
II we isolated the contributions of the rest-frame UV (~ 1220 
A), optical (~ 6570 A), Lya (1216 A), and Ua (6563 A). In this 
work, we measure the morphological parameters of these con¬ 
tributions. The LARS galaxies are hereafter referred to as Ln, 
where n ranges from 01 to 14 (see Paper 0). 


2.1. Morphological parameter estimation 

With the aim of quantifying the morphology of LARS galaxies, 
we calculated their sizes and performed non-parametric mea¬ 
surements of morphological parameters (see Appendix A and 
Pig. A.l for details). 

We calculated sizes, in terms of Petrosian semi-major axis 
(rP2 0 e.g.,[Lotz et al.|2004[[Lisker|2008[ ), circular Petrosian ra¬ 
dius ( [Petrokan[[1976| ), and~radii containing 20%, 50%, 80% of 
the flux (r 2 o, Go, I'so)- A comparison between these radii gives an 
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idea of the distribution of the light in the galaxy. We also esti¬ 
mated asymmetry (A), concentration (C), dumpiness (S), Gini 
coefficient (G), and second-order moment of the brightest 20% 
of the galaxy’s flux (M20, see e.g.,[Conselice||2003t|Lotz et al. 


|2004t[Striata et al.|2QQ7|[Micheva et al.|2013|). 

The asymmetry quantifies the symmetry of a galaxy with re¬ 
spect to a 180-degree rotation; the concentration describes how 
much the light is concentrated in the centre of a galaxy; the 
dumpiness measures the amount of small-scale structures within 
a galaxy; the Gini coefficient provides the information on how 
uniform is the fight distribution; M20 traces the spatial distribu¬ 
tion of any bright knots, and also off-centre clumps, its definition 
is very similar to that of C, but M20 is more sensitive to merger 
structures, such as off-centre components. 

We first ran the Source Extractor (SExtractor) software 
( Bertin & Arnouts 1996| ). It provided the galaxy centroid and 
the elliptical aperture, containing the entire galaxy and charac¬ 
terized by semi-major axis (sma) equal to rP20. The photometry 
was performed within this SExtractor detection aperture. 

We adopted configuration parameters like in [Bond et al.| 
( |200^ (DETECT_THRESH =1.65, DETECT_MINAREA=30, 
DEBLEND_MINCONT=l). They were optimized to provide 
morphological measurements in deep HST rest-frame UV ob¬ 
servations at z > 2. To prevent SExtractor from breaking up the 
clumpy, resolved z ~ 0 LARS galaxies into smaller fragments, 
we assumed a larger value of DETECT_MINAREA. This param¬ 
eter sets the number of contiguous pixels required for a detection 
to be accepted by SExtractor. We measured fluxes at SExtractor 
centroid within elliptical apertures, by using the ELLIPSE task 
in iraf.stsdas.isophote and within circular apertures, by using 
the PHOT task in iraf.digiphot.apphot. ELLIPSE and PHOT 
outputs served to infer sizes. A, and C at minimum asymmetry 
(CminA^^ as explained i n Appendix A and previously adopted in 
[Bershady et al. | ( [2000] ); |ConWice| ( [2003] ); |Micheva et al.|p013| ). 

The non-parametric measurements and signal-to-noise es¬ 
timations were performed counting the flux of pixels belong¬ 
ing to a segmentation map. We defined the segmentation map 
in two ways, one is an ellipse with semi-major axis equal to 
rP20 ( [Scarlata et al.|2007] ) and orientation given by SExtractor; 
one contains the pixels with surface brightnes s larger than the 
value at the Petrosian radius ( |Lotz et al.|2004| ) measured in the 
smoothed image (smoothed by a kernel of width rP20/5). We 
calculated M20, S and Gini coefficient by considering the pix¬ 
els within these segmentation maps. The Gini coefficients mea¬ 
sured in these two segmentation maps ar e denoted by G^ ^^Q and 


QSB-rp 2 {)s respectively. As described in Scarlata et al. ( 2007 ), 
QrP 2 Q defined to be consistent for redshift comparisons, thus 
we prefer it over throughout the paper when we com¬ 

pare with high redshift. 

To test our code, we applied it to template galaxies with 
known profiles and com pared the output to the results by 
[Bershady et al.| ( [2000| ) and jLotz et al.| ( |200^ . We recovered the 
expected values as described in Appendix A. 


2.2. Combination of morphological parameters 


As shown in jConselic^ ( [2003 1 ) and jLotz et al.| ( |20()4| ), combina¬ 
tions of morphological parameters can give information about 
galaxy history (e.g., star-formation and merging episodes). Eirst 
of all the rest-frame UV morphology is sensitive to the current 
star formation; the rest-frame opt ical traces the structure of the 
entire galaxy stellar population ( [Lee et aL 2Q13t Bond et al. 
2014| ). The combinations of parameters (see previous section) 


we adopted are. 


- Asymmetry vs concentration, asymmetry vs dumpiness, to¬ 
gether with dumpiness vs concentration, as presented by 
|Conselice| ( |2003] ) 

- Gini coefficient vs M20 bright-pixel moment, as presented in 
[Lotz et al.| ( |2Q()4| ). 

The concentration depends on the galaxy star-formation his¬ 
tory in the sense that a rapid gravitational collapse can produce 
high concentration. The presence of disk and intergalactic gas 
which cools onto the disk tends to produce a lower concentration 
value. Disk galaxies are character ized by 3 < C < 4, ellipticals 
by C > 4 ( [Bershady et al.|[2000| ). The asymmetry is sensitive 
to any feature that produces asymmetric fight distributions (e.g., 
star-formation knots, interactions, and mergers). It is commonly 
assumed also at h igh z that large asymmetry (A >0.38) in dicates 
a major merger ( [Aguirre et al.||2Ql^ [Conselic^|2003[ ). Spiral 
galaxies and systems composed of more than one component 
are characterized by A > 0.1. The dumpiness is sensitive to the 
presence of star-forming clumps as well, but background noise 
can make it difficult to detect low surface brightness regions and 
increase the appearance of the galaxy as a mix of clumps. The 
Gini coefficient can be strongly correlated with C. By definition, 
G = 1 means that the light is all concentrated in one pixel, G = 0 
that the light is equally distributed across the galactic body. In 
the case of a shallow fight profile, both G and C are low. When 
more than one clump contains a significant fraction of light, G 
can be much larger than zero, but C still low. M20 traces the 
spatial distribution of off-centre bright regions. 

In general, starburst and irregular galaxies are expected to 
have large A, large S, and intermediate C, merging systems and 
perturbed disks show large M20 and intermediate G. 


2.3. High-redshift simulation 

We simulated the observations of LARS galaxies (all at z < 0.2) 
at higher redshift by transforming their original science- and 
weight- map images (Paper I) according to the fo llowing steps 
(see also jOverzier et al.|2008t[Adamo et al.|2013[ ). 


1. The images were resampled preserving the flux (IDL frebin 
function). The size of the output image was defined by fixing 
the physical size of the galaxies. We chose mainly a z ~ 2 
sampling to be able to compare with the interesting results 
obtained by surveys of Lyman alpha emitters in the last re¬ 
cent 5 years ([Nilsson et al. 12009 [ [Guaita et al.|201Q[ [Hayes 


et al.|2010}[Nakajima et al.|2012[ Sandberg in prep.). Also, 

the size changes a little with redshift. _ 

2. Continuum subtraction ( [Hayes et'~ar |2009| ) was applied to 
the resampled images to generate rest-frame UV continuum 
and Lyof line, rest-frame optical continuum and Ha line im¬ 
ages. The fine images are in units of flux (erg sec“^cm“^), 
while the continuum images are in units of flux densities (erg 
sec“^cm“^A“^). 

3. The image pixel values were scaled based on luminosity 
distance and surface brightness di mming (i.e., Hubble & 
[Tolman|1935t[Bouwens et al.|2004[ ). 

4. Gaussian noise, corresponding to a certain simulated survey 
depth, was added to the resampled and rescaled images by 
running the MKNOISE task in iraf.artdata. To calculate un¬ 
certainty on galaxy sizes and morphological parameters we 
performed Monte Carlo simulations by repeating 100 real¬ 
izations of a noisy image. The noise applied was defined as 
the lOcr detection within a ~50 pixel (equivalent to a square 
aperture of -0.2” on a side for HST ACS optical Alters) area. 
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similar to the limits given for the HUDF (Hubble Ultra Deep 
Field, [Beckwith et al.|2006| ). We do not show simulations in 
which we only resampled the pixel scale to that of a ground- 
based telescope and instrument, because the main effect on 
continuum and line images was produced by survey depth 
and ground-based point spread function (PSF, see Sec.[^. 


To choose reasonable ranges of detection limits (Table to 
apply, we referred to the MUSYC ( MUlti-wavelength Survey 
by Yale-Ch ile) NB3727 narrow band ( jGuaita et al.|2Q10 [Bondj 


jet al.|2012] ), to the t riple narrow band byjNakajima et al. 


2012 |), 


CANDELS/HUDF ( jMcLure et al.|2013} bro ad band, and to the 
dual narrow-band survey by [Lee et aL ^012| ). 


3. LARS galaxies at z ~ 0 

To be able to compare the Lya, Ha, and continuum properties 
of LARS galaxies with those of high-z hya emitters (LAEs), we 
focused on the twelve LARS galaxies with EW(Lya) > 1 A 
as measured in Paper 11. Thus, we excluded from this analysis 
the two galaxies of the sample (L04 and L06) characterized by 
null Lya maps. LARS galaxies with integrated EW(LyQf) > 20 
A composed the subsample of LARS-LAEs (consisting of six 
galaxies). Various physical characteristics of the LARS galaxies 
(including their coordinates) are discussed in Paper 11. 

In Fig. [J] we present the RGB images of the twelve LARS 
galaxies. Most of the galaxies show localized knots of star for¬ 
mation superposed on extended rest-frame optical emission; the 
Lya emission is extended on larger angular scales (the Lya 
haloes). 

We investigated the properties of LARS galaxies in the con¬ 
text of other galaxy populations, to assess the fairness of our 
comparison. In Figs. and we show the location of the 
LARS galaxies in the half-light radius vs stellar mass (rso vs 
log(M=,/Mo)) and the half-light radius vs UV absolute magni¬ 
tude (r 5 o vs Mf/v) diagrams, to understand if the LARS galax¬ 
ies harbour stellar masses and UV magnitudes comparable to 
values in the literature. These diagrams have been designed for 
local galaxies, for which sizes could be easily measured in the 
rest-frame optical bands ( jShen et al.|2QQ3| ). However, measure¬ 
ments in the rest-frame UV could also be performed at high red- 
shift. Following the method described in Sec. we estimated 
the half-light radius as in the rest-frame optical and also in 
the rest-frame UV images. 

The high-z studies we adopted for comparison all performed 
size and morphological measurements by using HST images. 
These include. 


- Continuum-selected Lyman break galaxies (LBGs) at z ~ 3, 
with and without hya in emission ( Pentericci et al.|2010|); at 
z ~ 1, 2, an d 3 ( jMosleh et al.|2Qll|)7 at 1.5 z 3.6 ( |Law| 
et al.|2012|); at z ~ 1 .8 ( jLotz et al. 2004 ); z-drop outs at z ~ 7 


( jGrazian et al.|2Q12] ); high signal-to-noise z- and Y-drop outs 
detected in the Hubble Ultra Deep Field, UDF12 (|Ono et al. 

[20TT 


jet al. 


- Compact star-forming galaxies (cSFGs) at 2 < z < 3 (Barro 
. |2Q13| ). These authors have pointed out that, based on 


their number densities, masses, sizes, and star formation 
rates, z ~ 2 - 3 compact, star-forming galaxies were likely 
progenitors of compact, quiescent, massive galaxies at z < 2 


^ We to ok a dvantage of this codification of the Lupton 
jet ^ (|2004|) prescription to produce RGB 

http://dept.astro.lsa.umich.edu/~msshin/science/code/ 

Python_fits-image/ 


images: 


- Star-forming galaxies selected based on their B -z and z-K 

colour (sBzK, Yuma et al. 2012[ Lee et al. 2013); passive and 
star-forming galaxies se lected based on t heir B-z and z-K 
colour (pBzK and sBzK jLee et al.|2013] ) _ 

- Star-forming galaxies at z ~ 2 - 3 by [Law et aT ( |2Q12| ). 
These authors found a typical value of the Gini coefficient 
^QSB-rP 20 S _ Q 4 ^ Pqj. sources with the strongest Lya 
emission, characterized by ~ 1.5 x IQio MO 

- GOODS (Great Observatories Origins Deep Survey) and 

UDF (Ultra Deep Survey) z ~ 4 and GOODS z ~ 1.5 sources 
from the study of jLotz efak ( |2006|) __ 

- Sub-millimeter galaxies (SMGs, [Aguirre et al.|2013| ) 

- Narrow-band selected Lyman alpha emitters at z - 2.07 and 
z ^ 3.1 ( [Bond et al.[[20^[MT2l ) belonging to the MUSYC 
survey. We considered the stack of the z ~ 2.07 entire sample 
and of subsamples separated by photometric properties, UV- 
faint(UV-bright) with/? > 25.5(< 25.5), IRAC-faint(IRAC- 
bright) with f 3 . 6 //m < 0.57(> 0.57) /iJ, low-(high-)EW with 
EW( LyQf) < 66(> 66) A , red-(blue-)LAE with 5-/? > 0.5(< 
0.5) ( [Guaita et al.|201 ij ). 

- Narrow -band selected Ly man alpha emitters at z ~ 5.7,6.5, 
and 7.0 ( [Jiang et al.[2013| ), the first very-high-redshift sample 
where non-parametric morphological measurements were 
performed. 


The local-Universe studies, we adopted for comparison, in¬ 
clude. 


- Sloan Digital Sky Survey (SDSS) early- and late-type galaxy 
relations obtained from the analysis of images in the z band 
( [Shen et al.|2003] ) 


- Lyman break analogues (LBAs) at z - 0.2. These are 
local starbursts that share typical characteristics of high-z 
LBGs, such as stellar mass, metallicity, dust extinction, star- 
formation rate, an d physical size. We consider ed a sample 
of 30 LB As from [Qverzier et aT ( [2009 [ 2010 | ). They were 
characterized by a median absolute UV magnitude of -20.3, 
almost one magnitude fainter than typical LBGs. 


We found that the LARS galaxies occupy a quite wide range 
of r 5 o. Their rest-frame UV and optical sizes (Table are 
broadly consistent with LB As, LBGs, and SMGs. Their stellar 
mass tend to be larger than LAEs, consistent with LBAs and 
LBGs. However, there is an overlap in stellar mass between 
LARS-LAEs (M:, < 10^^ MG) and the most massive LAEs in 
the sample of [Bond et al.[ ( |2012| ). Also, LARS galaxies are less 
massive than cSFGs. The largest half-light radii characterize the 
LA RS ga laxies with the most distorted morphology (see also 
Fig. [B.2[ ). LARS M^/y magnitudes (and so star-formation rate, 
SFRf/v) are comparable with those of z - 2.07 LAEs and z > 7 
LBGs. There is an overlap with z > 5 LAEs. However, the mea¬ 
surements of LARS sizes in the rest-frame UV are larger than 
those of z > 7 LBGs. 

Therefore, LARS galaxies could be considered as LBAs, 
with size, stellar mass, and star-formation rate similar to 2 < 
z < 3 star-forming galaxies. 


3.1. Continuum morphology of LARS galaxies atz- 0 

Following the method described in Sec.[^ we estimated the non- 
parametric measurements for the LARS galaxies (Table [^. 

Combinations of morphological parameters (see Sec. |_ 

can give information about a galaxy’s star-formation history. 


Lotz et al. (2004) proposed a criterion for separating perturbed 
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Table 1: lOcr detection limits applied to high-z simulated LARS images 


F(Lya) 

m(rest-frame UV) 

F(Hcr) 

m(rest-frame optical) 

erg sec“^ 

AB 

erg sec“^ 

AB 

5E-19 

30 

2E-19 

29 

3E-18 

29 

6E-19 

28 

8E-18 

28 

IE-18 

27 

2E-17 

27 

3E-18 

26 

5E-17 

26 

lE-17 

25 


Notes. Noise corresponding to the detection limits given in the table was added to LARS continuum and line images. As LARS Lycr and 
Hcr(rest-UV a nd optical continua) images are in units of flux(flux density), the image depths are given in units of erg sec"^cm"^(AB magnitudes). 
The MUSYC ( Guaita et al.|2010 Bond et ari|2012| ) NB3727 survey implied a lOcr detection limit of F(LyQf)=5E-17 erg sec"^cm"^. MUSYC 
U,B(HUDF V606) lOcr detection limits were about 26(29.5). Assuming a NB3727 filter width and transmission profile, a source with m^ 5 (rest- 
frame UV)=30 and EW(Lycr) = 20 A is chara cterized by E(Lycr ) =5E-19 erg sec~^ cm~^. m^^(rest-frame UV) = 30,29,28 are consistent with 
HUDE09, GOODS, and GEMS survey depths ( [Bond et al.|2009| ). |Lee et al. \2012\ survey was characterized by a lOcr detection limit of 22.9 
in NB210 and 23.7 in K. CANDELS wide(deep) F160W lOcr detection limit was 25.8(26.5). Assuming a NB210 filter width, a source with 
m^fi(rest-frame optical) = 27 and EW(HQf) = 20 A is characterized by E(Hcr)=2E-19 erg sec"^cm"^. 





Eig. 1: Ealse-colour images of the LARS galaxies analysed in this paper. Red encodes rest-frame optical continuum, green rest-frame UV contin¬ 
uum, and blue shows continuum-subtracted Lycr. Scales in kpc are given on the side. Intensity cut levels are set to show details. 


disks or merging systems from normal galaxies, by studying 
local ultra luminous infrared galaxies (ULIRGs). The crite¬ 
rion identifies a region in the G vs M20 diagram, which is 
QSB-rP20s ^ _0.115 X M20 -f 0.384. For z < 1.2 galaxies 
obse rved in a rest-fram e optical band (4000 A) at HST resolu¬ 
tion, |Lotz et aL] (|2008|) proposed a slightly different relation to 


identify merging systems, > -0.14 x M20 -f 0.33. 

Also, |Conselice| ( [2003] ) distinguished the region where irregular 
or starburst galaxies were located in the and A-S planes. 

As seen in Fig. and Table ^ LARS galaxies, in particular 
the LARS-LAEs, tend to avoid the location of normal galax¬ 
ies and to occupy the region of perturbed disks or merging sys¬ 
tems and of irregular or starburst galaxies. The values of G, M20, 
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L. Guaita: LARS at low and high z 
Table 2: Size of the original LARS galaxies at z ~ 0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

LARS 

rP20"" 

rP20"'''" 


r20 

Lo 

TSO 


kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

rest-UV 







LOl 

2.63 

1.34 

0.81 

0.86 

1.20 

2.11 

L02 

2.42 

2.28 

1.49 

0.69 

1.28 

2.17 

L03 

1.40 

0.97 

0.77 

0.33 

0.63 

1.17 

L05 

1.71 

1.46 

0.84 

0.35 

0.69 

1.45 

L07 

1.13 

0.90 

0.81 

0.33 

0.60 

0.95 

LOS 

5.98 

3.35 

2.08 

1.67 

3.01 

4.96 

L09 

18.72 

1.11 

0.59 

2.73 

8.52 

13.81 

LIO 

4.73 

3.71 

3.09 

0.77 

2.08 

3.93 

Lll 

23.77 

18.56 

18.56 

6.12 

12.21 

19.41 

L12 

1.84 

1.20 

1.05 

0.41 

0.75 

1.76 

L13 

4.42 

3.70 

1.23 

1.95 

2.47 

3.75 

L14 

1.83 

1.71 

0.98 

0.67 

0.98 

1.52 

rest-optical 







LOl 

4.87 

2.15 

2.44 

1.28 

2.11 

4.01 

L02 

7.12 

2.23 

2.51 

1.08 

2.63 

5.62 

L03 

5.73 

1.98 

1.83 

0.72 

2.00 

4.59 

L05 

2.32 

1.79 

1.84 

0.41 

0.98 

1.79 

L07 

3.95 

2.89 

1.96 

0.62 

1.51 

3.18 

LOS 

5.28 

3.17 

3.29 

1.48 

2.63 

4.22 

L09 

20.54 

3.06 

2.07 

5.98 

10.27 

14.74 

LIO 

7.14 

3.93 

3.93 

1.08 

2.83 

5.57 

Lll 

19.13 

14.52 

14.52 

4.83 

9.82 

15.43 

L12 

3.68 

3.15 

1.20 

0.64 

1.46 

2.97 

L13 

11.41 

5.04 

5.35 

2.31 

5.30 

9.15 

L14 

3.72 

3.78 

1.22 

0.73 

1.28 

2.80 


Notes. Size measurements from elliptical and circular aperture photometry of the original LARS images. (1) LARS id, (2) Petrosian semi-major 
axis, (3) circular Petrosian radius, (4) Petrosian radius at minimum of asymmetry, (5) radius containing 20%, (6) 50%, and (7) 80% of the total 
flux. The measurements were performed in the bands corresponding to the rest-frame UV (either FI40 or FI50) and to the rest-frame optical (either 
F775, F814, or F850). The step in semi-major axis is 1 pixel (-0.02 kpc at z ~ 0.03) and in circular radius is 2 pixels (-0.05 kpc at z ~ 0.03). 


Table 3: Morphological parameters of the original LARS galaxies at z ~ 0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

LARS 

qvFIO 

QSB-rFmS 

M20 

ell 

Qcirc 

C“‘ 

QminA 

pixel 

A 

S 

rest-UV 











LOl 

0.70 

0.68 

-0.87 

0.61 

1.34 

1.94 

1.73 

251.40 

0.42 

0.16 

L02 

0.73 

0.72 

-0.79 

0.53 

3.14 

2.50 

3.49 

52.26 

0.50 

0.26 

L03 

0.63 

0.64 

-1.17 

0.30 

2.20 

2.71 

2.76 

83.18 

0.21 

0.11 

L05 

0.72 

0.68 

-1.26 

0.34 

3.11 

3.07 

3.27 

408.69 

0.33 

0.15 

L07 

0.62 

0.62 

-1.02 

0.40 

1.99 

2.28 

1.24 

696.00 

0.26 

0.09 

LOS 

0.62 

0.65 

-0.89 

0.33 

1.95 

2.36 

1.43 

21.98 

0.37 

0.29 

L09 

0.71 

0.75 

-2.33 

0.80 

1.76 

3.52 

3.35 

23.86 

0.21 

0.24 

LIO 

0.67 

0.68 

-1.52 

0.52 

3.53 

3.53 

2.76 

32.83 

0.33 

0.24 

Lll 

0.58 

0.79 

-1.22 

0.87 

3.40 

2.51 

3.40 

44.05 

0.30 

0.31 

L12 

0.67 

0.68 

-1.50 

0.34 

2.56 

3.15 

2.82 

583.81 

0.28 

0.04 

L13 

0.72 

0.67 

-0.77 

0.21 

1.23 

1.42 

2.61 

89.17 

0.33 

0.23 

L14 

0.74 

0.68 

-0.71 

0.27 

1.51 

1.78 

2.72 

622.68 

0.21 

0.02 

rest-optical 











LOl 

0.59 

0.61 

-1.09 

0.64 

1.72 

2.48 

2.53 

21.91 

0.29 

0.09 

L02 

0.59 

0.64 

-1.36 

0.69 

2.96 

3.58 

2.89 

4.25 

0.14 

0.31 

L03 

0.61 

0.62 

-2.16 

0.45 

3.32 

4.02 

3.42 

50.72 

0.09 

0.05 

L05 

0.64 

0.67 

-1.57 

0.73 

3.27 

3.22 

4.77 

64.96 

0.33 

0.15 

L07 

0.60 

0.61 

-1.68 

0.60 

3.28 

3.56 

2.97 

14.60 

0.24 

0.10 

LOS 

0.54 

0.55 

-1.05 

0.23 

2.10 

2.28 

1.78 

36.82 

0.31 

0.16 

L09 

0.62 

0.74 

-1.13 

0.83 

2.03 

1.96 

4.35 

11.38 

0.45 

0.34 

LIO 

0.57 

0.58 

-1.87 

0.54 

3.08 

3.56 

3.31 

9.95 

0.13 

0.15 

Lll 

0.50 

0.60 

-1.25 

0.85 

3.46 

2.52 

3.46 

18.70 

0.20 

0.23 

L12 

0.64 

0.63 

-1.90 

0.46 

3.70 

3.34 

3.35 

54.74 

0.33 

0.13 

L13 

0.58 

0.62 

-1.52 

0.55 

2.43 

2.99 

3.67 

3.98 

0.18 

0.30 

L14 

0.76 

0.72 

-1.32 

0.10 

2.92 

2.92 

3.01 

35.24 

0.32 

0.13 


Notes. Morphological parameters estimated for the original LARS galaxies, following the equations in Appendix A. (1) LARS id, (2) Gini 
coefflcient estimated within the flxed-size segmentation map, (3) Gini coefflcient estimated within the segmentation map built from the pixels with 
surface brightness larger than the value corresponding to that at the Petrosian radius, (4) M20, (5) SExtractor ellipticity 1-B/A, (6) concentration 
from circular apertures, (7) concentration from elliptical apertures, (8) concentration corresponding to the minimum of asymmetry, (9) signal-to- 
noise per pixel, (10) asymmetry, (11) dumpiness. As discussed in Appendix A, we could expect an uncertainty of - 10% in G, of < 5% in M20, 
and of - 20% in C. We could expect an even larger uncertainty in S. The difference in and is generally marginal. 
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Fig. 2: Half-light radius measured using elliptical apertures in the rest-frame UV {upper panel) and optical {lower panel) as a fu nction of stellar 
mass. LARS values at z ~ 0 from this work (red diamonds) are shown together with previous rest-frame UV measurem ents: [Overzier et aT] 
( |2009| ) individual LBA values estimated in the HST ACS/SB C F150LP and AC S/HRC F330W filter (open blue squares), [Yunia et a~ ( |2012| ) 
median estimations in the F850 W filter for sBzK (gre en star), [Bond et al.| ( [2012| ) LAE sample stack and subsamples (magenta squar e, blue and 
magenta sta rs, small black circle), [Aguirre et ar] ( |2013| ) individual SMG values measured in the HST FI lOW filter (big blac k dots), andjPentericcij 
jet al.| ( |20T0| ) average values of LBGs with and without Lycr in emission (green dots); and rest-frame optical measurements: [Overzier et al.| ( |2009j ) 
individua l LBA values estim ated in the HST (Wide Field and Planetary Camera2) WFPC2 /F606W and ACS/W FC F850LP filter (open blue 
squares), [Yuma et ^ ( |2012| ) median estimations in the F160W filter for sBzK (green star), [Mosleh et al. ( |2011| ) median values of UV-bright 
sources (GALEX-LB Gs at z ~ 0.6 - 1.5, LBG at z ~ 2.5 - 3.5, and c ontinuum-selected star-forming galaxies at z ~ 1.5 - 2.5, blue diamonds), 
[Aguirre et al.[ ( |201^ individual SMG values measured in E160W filter, [Law et al. (|2012| ) mean value of all the sample of star-forming galaxies at 
1.5 < z < 3.6 estimated in the E160W filter (yellow triangles), and |Barro et al.| ( [2013[) values for com pact star-forming galaxies at 2 < z < 3 also 
calculated in the E160W filter (green triangles). We also show the curve derived by lShen et al.[ ( |2003] ) for local SDSS early- and late-type galaxies. 
As these curves were obtained in z-bands for local galaxies, it is more meaningful to compare them to the radii in the rest-frame optical. However, 
for reference, we show them in the upper panel as well. The stellar masses are all corrected to Salpeter-IME values and the size measurements are 
all scaled to be comparable to half-light radii. LARS stellar masses were calculated in Paper II. L09 and LI 1 are outside the graph, due to their 
half-light radius larger than 6 kpc. 


C, and A, we calculated for LARS galaxies, are consistent with ter, characterized by one of the smallest stellar masses and the 

the ones measured by [Overzier et al.[ ( [201Q[ ) for LBAs. Even if largest among the LARS-LAEs. 

our sample is just composed of twelve sources, we do not see 

any significant dependency between G and EW(LyQf). LOS is 

the most massive of the LARS-LAEs, but equally concentrated 

within the segmentation map. L02 is the largest-EW(LyQf) emit- 
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L. Guaita: LARS at low and high z 
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Fig. 3: Half-light radius measured using elliptical apertures in the rest-frame UV image as a function of the absolute rest-fram e UV magnitude. 
LARS measurements from this work (red diamonds) are shown together with the literature estimatio ns by |Overzier et al.|f2010| ) corresponding to 
the medi an value of their z < 0.3 LB As, observed in the near-infrared bands, by | Jiang et al.|p013] ) for a sample of z - 5.7,6.5,7.0 LAE s (cyan 
dots), by [Bond et al.| { |201^ for the stack and subsamples of LAEs at z - 2 - 3 (black c ircles, magenta sq uares), by |Grazian et al.|p012| ) which 
measured SExtractor half-light radii for a sample of z-drop outs (green triangles), and by |Ono et al.| ( |2013] ) for a sample of high signal-to-noise z- 
and Y-drop outs detected in Hubble Ultra Deep Eield, UDE12 (black stars). 


3.2. Lya morphology of LARS galaxies atz- 0 

One of the goals of our work was to quantify and compare the 
morphologies of LARS galaxies in Lya and in the continuum. 
We present morphological parameters measured in the continua 
and in Lya of LARS images in Figs. and [6 M20, 

concentration, and ellipticity are smaller; while dumpiness and 
asymmetry are generally larger in Lya than in the rest-frame 
UV continuum. The LARS-LAEs tend to be characterized by 
the highest concentration, lowest asymmetry, and lowest dumpi¬ 
ness in Lya. and M20 measured in the rest-frame optical 
are consistent with the values measured in Lya. 


4, LARS galaxies as seen at z ~ 2 


We applied the procedure described in Sec. |2.3| to simulate 
LARS galaxies at z ~ 2. We named the high-z simulated galax¬ 
ies as z2LARS and the subsample of Lya emitters as z2LARS- 
LAEs. We estimated sizes and calculated morphological param¬ 
eters (Tables E-E an d [9p^ in the same way we did for the 
original images in Sec. [2.1 


The purpose of this test was to understand whether we could 
expect to detect LARS-type galaxies and LARS-type Lya haloes 
in current high-z surveys. In particular, we wanted to understand 
how galaxy size and morphological parameters changed when 
varying the survey depth (Table [^. The results show that, in a 
sufficiently deep survey, faint galaxy structures in between bright 
knots remain connected together and SExtractor is able to detect 
just one source (the entire galaxy) in the image. In a shallower 


survey, the faint connecting structures tend to be lost in the noise 
and a galaxy appears to be composed of separated clumps. In that 
case SExtractor identifies more than one source and photometry 
is performed by locatin g the pho tom etric aperture around the 
brightest clump. In Eigs. B.1||B.2[ and B.3 we show how LARS 
galaxies would appear if detected in the deepest continuum and 
line surveys simulated here, while Eigs. |B.4 |B.5| and B.6 show 
the results for shallower surveys. In Appendix |C| we present the 
corresponding surface brightness profiles. 


In the following sub-sections, we describe the detection of 
LARS galaxies in the simulated surveys with lOcr detection lim¬ 
its presented in Table In the first sub-section, we give details 
on the detection of LOl as an example. We proceed to describe 
the cases of the LARS-LAEs and of the galaxies with the faintest 
Lya emission. Then, we explain the variations in size and el¬ 
lipticity versus dumpiness owing pixel resampling and survey 
depth. In Sec. |4.6[ we quantify the morphology of z2LARS and 
compare with high-z observations from the literature. 


4.1. Detection of LOl in high-redshift surveys 

In Eig.[7] we show the rest-frame UV, Lya, optical, and Lia im¬ 
ages of LOl simulated to be at z ~ 2 (z2L01) as they would 
be observed in the deepest surveys probed here; SExtractor de¬ 
tection apertures are over-plotted. The detection parameters we 
adopted are sensitive enough that (at the deepest simulated sur¬ 
veys) this galaxy is detected as a single source. As described in 
detail in Paper I and II, LOl consists of a bright UV star-forming 























































L. Guaita: LARS at low and high z 
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Fig. 4: Combinations of rest-frame optical morphological measurements, used in the literature as diagnostics of galaxy past and current history. 
QSB-rP 20 s yg jy[2o {upper left), A vs {upper right), and A vs S {lower right). The lower /g ft panel shows stell ar mass vs the rest-frame 


UV rPiQs ^ typical value of G^^ 


0.4 for the strongest Lya emitters of the sample by 


Law et al. 


<2012 1 is reported as a red star. 


The twelve LARS galaxies analysed here are presented as squares, LARS -LAEs are rounded by open circles. The colour scale corresponds to 
EW(LyQf). Eor co mparison, green sta rs correspond to the|Frei et al.|(|1996| ) sample of normal galaxies and light blue triangles correspond to the 
ULIRG sample of|Borne et al.|(|2000|) as processed by |Lotz et al.H2()04| ). Black triangles correspond to a sample of starburs t galaxies presente d in 


Conselice 

2003 

Lotz et al. 

2008 


I respectively (see text). 


centre with an extended tail, also seen in Her and in the rest- 
frame optical. The Lyer emission is coincident with the bright 
UV knot and extends in a fan-like structure possibly indicating 
the presence of an expanding bubble. The main features of emis¬ 
sion (dark red pixels in Fig. and absorption (white pixels), 
observed in Lycr thanks to the HST resolution and the careful 
continuum subtraction presented in Paper II, are clearly visible 
in the z ~ 2 simulation as well. However, the extremely detailed 
Lycr structures close to the centre of the galaxy (see Paper I Fig. 
1) are not visible. The last panel of Fig. |7] shows LOl Lycr image, 
convolved with a ground-based seeing. From the ground LOl 
Lyof morphology would appear smoothed. 

We show the surface brightness profiles of z2L01 in Fig.[^ 
The rest-frame UV and Lycr profiles (left column panels) are pre¬ 
served when observed in a survey with sensitivity deeper than 


"^^rTst-uv ~ F(LyQf)^^"^ = 8E-18 erg sec“^ cm"^. However, 

on scales larger than 4 kpc, the profiles are indistinguishable 
from the background noise. In shallower surveys, the profiles 
start to be affected by the simulated-survey noise on smaller 
scales and z2L01 could not be detected by adopting a SExtractor 
detection threshold, DETECT_THRESH = 1.65. Therefore, size 
and morphological parameter measurements could not be per¬ 
formed either. We define E(LyQf)^^"^ as the limits 

for detection and morphological parameter measurement. These 
limits are ^rZt-opticai ~ E(Hq')^*"^ = 3E-18 erg sec“^ cm“^ 

for LOl rest-frame optical and Her. 

The lower left panel of Eig.[8 shows that the rest-frame UV 
continuum profile is steeper than the Lyer profile. The lower right 
panel shows that the rest-frame optical continuum tends to be 
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L. Guaita: LARS at low and high z 



S(rest-UV} l-B/A(rest-UV) 


Fig. 5: Non-parametric measurements performed in the Lya images versus the ones performed in the rest-frame UV. From the upper left to 
the lower right: M20, A, S, and SExtractor ellipticity (1-B/A, where A and B are the semi-major and semi-minor axes of the 

detection ellipse). The dashed line indicates the 1:1 relation. The numbers reported in each panel correspond to the Spearman test coefficient, r, 
and probability, p, of uncorrelated datasets, r = 0 indicates no correlation, r = 1(-1) indicates direct(indirect) proportionality. 


shallower than the rest-frame UV and the Her, and more similar 
to the Lya profiles. 

4.2. Detection of the LARS-LAEs in high-redshift surveys 

The LARS-LAE galaxies show more than one bright knot, con¬ 
nected by filaments, in the continua. They also show an intense 
Lya emission close to their centres and Lya structures in their 


outskirts. The Lya emission is accompanied by regions of ab¬ 
sorption. As a typical trend, in increasingly shallower surveys 
the filaments, seen in the continua, show lower surface bright¬ 
ness, while their Lya emission become increasingly localized in 
the galaxy centre. Only LI4, the galaxy brightest in Lycr, could 
be detected in the shallowest Lya survey probed here. The mag¬ 
nitude and flux limits for z2LARS-LAE detection are shown in 
Tablell 
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S( rest-o pti cal) 1-B/A( rest-optica I} 


Fig. 6: Non-parametric measurements performed in the Lya images versus the ones performed in the rest-frame optical. From the upper left to 
the lower right we show the same parameters as in Fig.j^ 


4.3. Detection of the Lya-faint z2LARS galaxies 

In general these LARS galaxies are bright enough in UV, opti¬ 
cal, and Ha (see Paper II) to be detected as a sing le source in our 
SExtractor run. The exception is LI3 (Fig. 1^: this galaxy is 
detected as two possibly blended sources (dashed-line aperture 
in Fig. in UV and Ha, and photometry was performed lo¬ 
cating the aperture on th e (brightes t) ri ght c lump. In the shallow 
simulated surveys (Figs. B.4| B.5 and |B.6| ) low surface bright¬ 
ness filaments connecting the main continuum knots disappear 
into the background noise making the sources appear to be com¬ 


posed of several clumps. In increasingly shallow surveys, fewer 
clumps could be detected. For example, only the brightest clump 
is detected for LI3 at mrest-uv = 29 and oirest-opticai = 26. Only 
L09, the brightest galaxy in Ha, could be detected in the shal¬ 
lowest Ha survey probed here. 

In general these galaxies could not be detected in Lya sur¬ 
veys shallower than F(Lya) = 3E-18 erg sec“^ cm“^. The galax¬ 
ies presenting the strongest Lya absorption (white pixels in Fig. 
B of the sample could be detected only in the deepest Lya 
survey probed here. 
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Fig. 7: Simulated z ~ 2 rest-frame UV, Lycr, rest-frame optical, and Her emission for LOl as it would be observed in the deepest surveys probed 
here (first four panels starting from the left). The last panel shows the original-pixel-scale Lycr image convolved with a Gaussian kernel, resem¬ 
bling a ground-based seeing Point Spread Function, PSF, of 1.32”. Every panel is 20x17 kpc wide. The reddish ellipses indicate SEx apertures, 
corresponding to the assumed detection parameters: DETECT .THRESH =1.65, DETECT _MINAREA=30, and DEB LEND _MINCONT=l from 
Bond et al. (2009). The log colour scaling is chosen to show a visually consistent background noise. 






Fig. 8: Normalized surface brightness profiles of z2L01. The semi-major axis of an elliptical aperture is indicated as sma. Black points with error 
bars correspond to the surface brightness profile of the original LARS images in the rest-frame UV, optical, and Lycr line. The red squares represent 
the profiles at the limits of detection. The profiles indicated with circles(triangles) correspond to the deepest(intermediate) depth surveys. In the 
upper left panel we show the original rest-frame UV profile and the one at a depth of nirest-uv - 30, 29(cyan) and 28(red). In the upper right 
panel the original rest-frame optical profile, the one at a depth of rurest-opticai - 29,28,27(cyan), and 26(red). In the lower left the original Lycr 
profile, the one at a depth of F{Lya) - 5E-19, 3E-18(cyan), and 8E-18(red), together with the deepest-survey rest-frame UV profile (green circle 
and line). In the lower right the Her profile at a depth of F{Ha) - 2E-19, 6E-19(cyan), and 1E-I8(red), together with the deepest-survey rest-frame 
optical profile (green circle and line) and the deepest-survey rest-frame UV profile (black line with error bars). 


4.4. Size As proposed in Paper 0, we define the quantity ^ to estimate 

the size of Lya with respect to the size of Her (eq. and contin¬ 
uum (eq. emission, 


The Petrosian radius measured within elliptical apertures 
(rP20^^^) is the quantity adopted for the comparison of Lycr ver¬ 
sus continuum size estimations (Appendix A and Paper 0). The 
rP20^^^ estimated in the z2LARS continuum and Lya images 
could vary by up to 20% in median (Fig. [^. The variation de¬ 
pends on the specific morphology. 


^(LyalHa) = 


rP20^^\Lya) 

rP20"^^(ifQf) 


( 1 ) 


^(LyalUV) = 


rP20^^^(Lya) 
rP20^^Krest - UV) 


( 2 ) 


A ^ larger than 1 indicates that Lya photons extend to larger 
scales than UV and Her due to neutral hydrogen scattering; the 
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L. Guaita: LARS at low and high z 
Table 4: Limits of detection for LARS galaxies simulated to z ~ 2 


ID 

r^^rest-UV 

rarest-optical 

F(Lya) 

E(Ha) 


AB 

AB 

erg sec“^ cm“^ 

erg sec“^ 

LOllae 

28 

26 

8E-18 

3E-18 

L021ae 

29 

27 

8E-18 

lE-18 

L03 

29 

25 

3E-18 

3E-18 

L051ae 

27 

26 

8E-18 

3E-18 

LOTlae 

28 

26 

2E-17 

3E-18 

LOSlae 

28 

25 

8E-18 

3E-18 

L09 

28 

27 

5E-19 

lE-17 

LIO 

28 

25 

5E-19 

lE-18 

Lll 

28 

26 

3E-18 

3E-18 

L12 

26 

26 

3E-18 

3E-18 

L13 

28 

27 

8E-18 

3E-18 

L141ae 

27 

27 

5E-17 

lE-17 


Notes. Almost all LARS galaxies would be detected in Lya at a lOcr detection limit depth of 3E-18 erg sec"^ cm"^. The exceptions are L09 
and LIO, which show strong absorption in their Lycr images. LARS-LAEs would be detected in shallower surveys. Also, the majority of LARS 
galaxies would be detected in Her at a depth of 3E-18 erg sec"^ cm"^. 



size in the z2LARS images size in the original LARS images 


Eig. 9: Left: rP20^^^ (kpc) measured in the original LARS images as a function of the value from the high-z simulated ones. The 1:1 relation 
is shown as a dashed line. Right: Asize = ^ function of rP20^^^ measured in the original LARS images. The black 

dashed line indicates A = 0; the dashed blue, green, and red lines the median A values for Lycr, rest-frame UV, and optical. Open circles indicate 
LARS-LAEs; blue, green, and red dots correspond to the measurements in Lya, rest-frame UV, and optical images. 



$ in the z2LARS images ^ in the z2LARS images 


Fig. 10: ^(Lycr/rest-UV) (left) and ^(rest-optical/rest-UV) {right) measured in the original as a function of the regridded LARS images. Open 
circles indicate LARS-LAEs. The dashed line shows a 1:1 relation. 
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L. Guaita: LARS at low and high z 


galaxy presents a Lya halo. The left panel of Fig. [^compares ^ 
measured in the high-z simulated image with that of the original 
LARS image. Even if the rP20^^^ estimated in the high-z simula¬ 
tion is larger than in the original image, the ^ values remain con¬ 
sistent within the errors. The ratio between ^(Lyof/UV) measured 
in the z2LARS and original LARS images is 1.1 + 0.8 on aver¬ 
age. The ratios ^ (Lya/UV) or ^(Lya/Ha) vary between 1 and 5 
among the LARS galaxies, implying that Lya haloes are com¬ 
mon in the LARS sample. However, the total extension of the 
haloes depends on the depth of the simulated survey and the lim¬ 
its of detection in Lya, UV, and Ha (Sec. |4.1| 4.2, and 4.3). The 
sizes of z2LARS for each simulated survey probed here are pre- 
sented i n Tables [ 7 ] and[^ For comparison |Finkelstein et al.| 
( 2011b| ) reported size measurements of three LAEs at z ~ 4.5 of 
1-2 kpc in the rest-frame UV and of 1-3 kpc in Lya. Similarly 
Bond et al. ( 2010| ) estimated 2 kpc both in UV and in Lya. 

In Fig. 10 we also show 


^(optical/UV) = 


rP20^'^\rest - optical) 
rPlO^^Krest - UV) 


( 3 ) 


From the right panel of the figure, we see that the ratio be¬ 
tween ^(optical/UV) measured in the z2LARS and in the original 
LARS images is 1.1 + 0.3 on average. 


4 . 5 . EHipticity and dumpiness 


As noticed in [Bond et ah] ( |2009| ), it is not possible to easily 
understand if clumps in high-z broad-band images are merg¬ 
ing components rather than star- forming regions connected by 
low surface brightness structures. [Gronwall et al.| pOl 1| ); Bond] 
|et al.| ( |2012| ) observed that the rest-frame UV emission from a 
typical LAE was neither smooth nor spheroidal; its ellipticity 
was about 0.6, suggesting the presence of elongated structures 
due to merging activity or clumps of star formation. However, 
[Shibuya et al. ( 2014[ ) found that the LAEs in their sample with 
EW(Lya) > 100 A tended to be characterized by a small ellip¬ 
ticity both in the rest-fram e UV and optical. This is s upported 
by the theor etical results of Verhamme et al. ( |2012| ) and |Laurs^ 
|et al.| ( [2009] ), which showed that Lya photons could more easily 


escape from face-on disks (ellipticity ~ 0), generally character¬ 
ized by low column density of HI along the line of sight. 

LARS galaxies do not show any clear trend between contin¬ 
uum ellipticity and Lya equivalent width (Table [^. After pixel 
resampling the rest-frame UV, ellipticity varied by up to 60% 
for some of the LARS galaxies and no trend remains between 
continuum ellipticity and Lya equivalent width for the z2LARS 
galaxies. 

In shallow surveys, bright knots of star formation, seen with¬ 
out low surface brightness connectors, could appear aligned 
and lead to the galaxy appearing more elongated. Eigs. pT]and 
T2 show ellipticity versus dumpiness for a few z2LARS for 
which these measurements are significant. They also show how 
those parameters change when measured in the range of surveys 
probed here. Quantitatively z2LARS galaxies tend to have lower 
ellipticity in Lya than in the rest-frame UV, optical continuum, 
and Ha. Also, z2LARS-LAEs tend to show lower S than the 
other LARS galaxies in Lya. 

On average the ellipticity values increase in a shallower and 
shallower survey when measured in the rest-frame optical. We 
evaluate the correlation between depth and ellipticity or dumpi¬ 
ness by calculating the difference between the morphological pa¬ 
rameter measured in shallow surveys and in the deepest one. In 
the sample of LARS galaxies, 70% show a Spearman probability 


p ~ 0.0 and coefficient r-1 for the correlation of ellipticity. 

Also, 85% of the LARS galaxies show a correlation between 
depth and S estimated in the rest-optical. Eor the z2LARS-LAEs 
the ellipticity increases when measured in the rest-frame UV as 
well. Therefore, some of the high-z observations of large ellip¬ 
ticity and clumpy systems could be explained in terms of survey 
depth. Only a small change in S and ellipticity is seen in the Ha 
images. In Lya we generally measure a significant increase of 
S in shallow surveys. Of the LARS galaxies, 85% are charac¬ 
terized by a correlation between depth and S when measured in 
Lya. Eor a few galaxies the Lya ellipticity also increases with 
the shallowness of the survey. 


4 . 6 . Morphological parameters of LARS atz-2 


As in the case of the original LARS images (Sec. |3.1| ), we 
performed non-parametric measurements of morphological pa¬ 
rameters for z2LARS. To compare with high-z literature obser¬ 
vations, we made used of the G vs M20 diagram. We chose 
qtpiq shows a better consistency for high-z comparisons 

( [Scarlata et al. 12007] ). As shown in Appendix A, the pixel resam¬ 
pling has little effect on the measurements of and M20 for 
LARS galaxies. The change of parameters owi ng resampling has 
been studied in the literature by some groups (|Lotz et al.|2006[ 


Qverzier et al.||2Q10}|Huertas-Company et al.|[2014[ [Petty et al.| 


2014). The changes for irregular galaxies are in agreement with 
our findings for LARS galaxies. The non-parametric measure¬ 
ments for all the simulated surveys probed here are presented in 
Tables and The left panel of Eig. shows that 

z2LARS have morphology consistent with that of merging sys¬ 
tem, l ike z ~ 2.5 SMGs and some of the LBGs studied by |Lotz 
( |2004| ) . They tend to be separated fro m the location of z 


et al. 


sBzK and pBzK galaxies ( |Lee et al.|2013| ). Some z2LARS-LAEs 
show a larger G^^^^ per M20 value than the other z2LARS. 


Lotz et al. ( 2006| ) built a criterion to identify merging sys¬ 


tems based on their rest-frame UV. Since the rest-frame UV 
tends to be more disturbed than the optical, they defined two 
main regions in the G-M20 diagram. The condition M20 > -1.1, 
typical of well-separated double or multiple bright nuclei, was 
used to identify major mergers, while the condition M20 < -1.6, 
G > 0.6 to identify bulge-dominated systems. The left panel of 
Eig. shows that z2LARS galaxies have morphologies simi¬ 
lar to that of high-z LAEs and LBGs in the rest-frame UV, and 
tend to be less compact than SMGs. In this diagram LOl, L08, 
and L14 have morphologies consistent with that of major-merger 
systems. The non-LAE z2LARS galaxies, similar to the major¬ 
ity of the star-forming galaxies at high redshift, are character¬ 
ized^^ intermediate M20 values. In the right panel of Eigs. 
andf^ we present G^^^^ vs M20 measured in the Lya images. 
The subsample of z2LARS-LAEs tend to have lower G^^^^ and 
lower M20 in Lya than in the rest-frame optical and UV. Some 
of the z2LARS non-LAEs tend to have a more distorted (com¬ 
posed of more than one clump) morphology in Lya than in the 
continua. It can be seen in Tables [9p2| that the morphology of 
LARS galaxies does not change significantly in the rest-frame 
optical when non-parametric measurements were performed in 
simulated surv^s with depths comparable with the limits of de¬ 
tection (Table 1^. Eor a few z2LARS, G^^^^ and M20 vary in a 
way that they approach the dashed line drawn in Eig. In the 
rest-frame UV, G^^^^ and M20 become characterized by larger 
uncertainties, while in Lya z2LARS become increasingly com¬ 
pact and less composed of multiple structures. 
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L. Guaita: LARS at low and high z 




Fig. 11: Ellipticity as defined by SEx (1-B/A, where A and B are the semi-major and minor axes of the detection ellipse) as a function of dumpiness 
for the rest-frame UV and optical continuum. Diamonds represent z2LARS measurements in the deepest survey, triangle and circles in increasingly 
shallower surveys. Dashed line is a 1:1 relation drawn to aid the eye. 




Eig. 12: Ellipticity versus dumpiness as in Eigj^for measurements performed in Lycr and Her images. 


5. Stacking of regridded LARS images 

As extensively described in Paper 0 and II, a significant fraction 
of the Lya photons in LARS galaxies is emitted from haloes. 
These haloes begin in the inner few kpc and extend outward to 
scales larger than those characterizing localized star-forming re¬ 
gions. Lya maps were shown in Paper II as well as maps of 
Lyof/Hof ratio and dust-reddening maps. Since Ha photons are 
emitted directly from the HII regions (i.e., they do not scat¬ 
ter), a value of the Lya [Ha ratio that exceeds that of case B 
recombination is most probably related to HI scattering. Dust 
reddening and HI scattering together can contribute to the situa¬ 
tion in which the number of Lya photons observed on a partic¬ 
ular sight-line is reduced. The re-processing of ionizing photons 
(e.g., [Humphrey et~ar]|2008| ) in a region of the galaxy different 
from that of the star formation could contribute to an additional 
production of Lya photons and a possible extended Lya emis¬ 
sion. In this case, however, Lia radiation will also be produced at 
large radii by the same recombinations that make Lya. The avail¬ 


ability of Lia images, together with those of Lya, has favoured 
an interpretation in which HI scattering plays a key role in Lya 
emission, but it has proven difficult to disentangle the contri¬ 
bution of other factors. However, HI scattering would also pro- 
duce a certain degree of polarization in th e Lya emission (e.g., 
[Humphrey et al. [20ISbj[Hayes et al.|2QlT] ). At high redshift, Lya 
haloes have been detected only in a minority of cases, indicating 
that their detection is challenging. 


5.1. Literature results 


[Rauch et al.[ ( [2008[ ) stacked Lya spectra of a sample of faint 
(V > 25.5 for 80% of the sample) spectroscopically detected 
galaxies at 2.67 < z < 3.75, reaching a Lya surface brightness 
limit of ~ 2E-19 erg sec“^ cm“^ arcsec“^ at Icr. The resulting 
Lya profile extended up to ~30 kpc from the centre. [Matsuda 
et al. ( 2012[ ) stacked narrow-band images of large subsamples of 


LBGs, separated based on their environment. They found that 
galaxies in overdense regions tend to show Lya emission on 
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Fig. 13: Left: vs M20 measured in the rest-frame optical images of z2LARS (big diamonds). Open circles i ndicate z2LARS- LAEs. Data from 

the literature are presented as s mall symbols: L BGs detected in the Hubble Deep Field-North and meas ured by|Lotz et al.|(|2004| ) (green squares), 
passive and star-forming BzK ( |Lee et aL|2013| grey diamonds and yellow squares), and z ~ 2.5 SMGs ( [Aguirre et al.|2013| black circles), where 
the measurement s are all performed on HST images. The dashed line indicates the separation between normal galaxies and merging systems by 


Lotz et al. 


2008| ). Right : G^^^^ vs M20 measured in the Lycr images of z2LARS (big diamonds). The colour coding for LARS galaxies is as used 


throughout this paper. The dashed red line represents the rest-frame optical (left panel) separating region, drawn to guide the eye. 
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Fig. 14: Left: G^^^^ vs M20 measured in the rest-frame UV images of z2LARS (big diamonds). Open circles indicate z2LARS-LAEs. Small 
symbols are data f rom literature: LAEs and LBGs at z > 6 ( jJiang et al.|2013 cyan diamonds), GOODS z ~ 1.5 and G OODS/UDF z - 4 gala xies 
( jLotz et al.|2006[ grey circles and yellow tri angles), mean values for z ~ 3 LBGs with and without Lycr emission ( jPentericci et al.|2010| red 
triangle and black cross), and z ~ 2.5 SMGs ( [Aguirre et al.|2013"] blac k circles), where t he measurements are all performed in HST images. Eor 
comparison we also show the local galaxies regridded to be at z ~ 1.5 by [Lotz et al.[ ( [2006] ): an elliptic al (open circle), an Sb (open triangle), an Scb 
(open square), and a merging galaxy (open star). The dashed lines indicate the criteria introduced by [Lotz et al.[ ( [2006] ) to identify major mergers 
and bulge-dominated systems. Right: G^^^^ vs M20 measured the Lycr images of z2LARS (big diamonds). The dashed green lines represent the 
rest-frame UV (left panel) separating regions, drawn to aid the eye. The colour coding for LARS galaxies is the same as in Eig. 
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scales up to 40-60 kpc, and are larger than those of isolated 
galaxies. This c ould imply that Lya haloes follow the dark mat¬ 
ter distribution. [Steidel et al.[ ( |MTT] ) (hereafter Sll) also stacked 
Lya images of subsamples of LBGs at z ~ 2.65, segregated by 
Lya equivalent width. The galaxies came from surveys of three 
different fields, which were also characterized by over-densities. 
The stack of the entire sample reached a Lya surface brightness 
limit of ~ IE-19 erg sec“^ cm“^ arcsec“^ (Icr), while the LAE 
(EW(Lya) > 20 A) subsample reached ~ 2.4E-19 erg sec“^ 
cm“^ arcsec"^. They found that the stacked Lya profiles had a 
characteristic size up to 9 times larger than the stellar continua. 


Their detected Lya emission extended up to 80 kpc from the 
centre of the stacked source. 


By stacking images of 187 narrow-band selected Lya emit¬ 
ters at z - 2.07, [Eeldmeier et al.[ ( [2013[ ) found no evidence of 
extended Lya emission, but did observe a 5-8 kpc halo in a 
stack of about 200 LAEs at z - 3.1. These field-galaxy stacks 
reached a Lya surface brightness Icr limit of just 9.9E-19 erg 

2 07 and 6.2E-19 erg s ec~^ cm“^ 
3.1. Recently, Momose et al. ( |2014) (hereafter 
Ml4) succeeded in detecting extended haloes, by stacking over 
3500, 300, and 350 narrow-band detected LAEs at z ~ 2.2, 


sec 


cm ^ arcsec ^ at z 


arcsec ^ at z 
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L. Guaita: LARS at low and high z 


~ 3.1, and z - 5.1 from Subaru surveys (not necessarily in over- 
dense regions). They reached a Lya surface brightness Icr limit 
of 1.6E-20, 1.7E-19, and 5.5E-20 erg sec“^ cm“^ arcsec"^, re¬ 
spectively. Both these last two studies performed L AE selection 
in deep (5cr detection limit of 25 AB for Eeldmeier; |Guaita et al.| 
|2010t[Gronwall et al.|2007| ) and very deep (5(T dete ction limit of 
25.1-25.7 AB for Momose; Nakajima et al.|20T2 ) narrow-band 
images. Also, they treated the sources of uncertainty in a very 
careful way. However, the detected haloes (r^ ~ 8 kpc at z = 2.2) 
were not as extended as the ones claimed by Sll (r^ ~ 25 
kpc at z = 2.65). A few attempts to detect Lya haloes from 
individual high-z LAEs were performed by |Bond et (|2010|) 
and|Einkelstein et ak] ( |2011a| ) in HST-filter images. |Bond et al.| 
( 2010| ) presented a morphological study of LAE Lya emission at 
z - 3.1 (E502N filter, Icr detection limit of 3E-17 erg sec“^cm“^ 
arcsec“^ for 1 arcsec source). Their images were very shallow 
and did not show significant extended haloes (^(LyajUV) ~ 1, 
Table [5]). An estimation at z ~ 4.5 was presented in |Einkelstein 


|et al.|^^lla| ) (E658N filter, Icr detection limit of about 2E-18 
erg sec“^cm ^ for 1 arcsec source). They studied Lya emission 
from three spectroscopically confirmed LAEs and found evi¬ 
dence of Lya haloes in two of them (^(LyajUV) ~ 1.5). 


5.2. Stacking procedure and stacked surface brightness of 
LARS at high redshift 


By following the steps listed in Sec. we simulated the ap¬ 
pearance of LARS galaxies at z ~ 2 and z - 5.7. After adding 
noise to the simulated high-z images, we stacked the observa¬ 
tions of individual galaxies the same way it is done at high red- 
shift to increase the signal-to-noise. We used the IMCOMBINE 
task in iraf .images.immatch to (average) stack every galaxy at 
the position of SExtractor centroid. We generated an average 
(LARSaverage), median stack of all the twelve high-z LARS 
galaxies (like in Mil), and also an average stack of the six 
LARS-LAEs. It was only meaningful to look at the stacked 
profiles up to 12 kpc, which was the physical scale probed by 
the ACS/SBC detector common to all LARS images. However, 
narrow-band observations at high redshift are able to probe much 
larger scales. As the median stack provided very similar results 
as the average one, we only reported results from the LARS- 
LAE and LARS average stacks in the tables and figures. We 
chose to combine individual galaxy frames from simulated sur¬ 
vey depth of E(Lya)=5E-19, 3E-18, 2E-17, and lE-15 erg sec“^ 
cm^. The central value of E(Lya)=2E-17 is comparable to re¬ 
cent ground-based narrow-band surveys. The sha llowest value 
was chosen to match the depth of Bond et al. ( |2010| ), in which we 
did not detect any Lya emission from the stacks. The final depths 


depended on the number of sources in the stack, ¥(Lya)l Vl2 for 
LARSaverage and ¥(Lya)|^^6 for LARS-LAEs. We applied the 
same code we ran on the individual z2LARS galaxies to obtain 
surface brightness profiles (Eig. [m and sizes (Table of the 
stacks at z ~ 2 and z - 5.7. The LARS-LAE Lya profile was 
the most peaked in the centre, and was also the most affected by 
background noise at large radii. The Sersic indices of the LARS- 
LAE stack profiles were found to be 3.4 and 2.1, while the ones 
of the LARSaverage stacks are 2.5 and 2.0 at z ~ 2 and z - 5.7, 
respectively. 


We find that a depth of E(Lya) = 5E-19 erg sec“^ cm“^ 
(~1.4E-19 erg sec“^ cm“^ after stacking twelve sources) enables 
us to recover Lya haloes in the stacks of LARS galaxies at both 
z ~ 2 and z - 5.7. A depth of E(Lya) = 2E-17 erg sec“^ 
cm“^ (~5.8E-18 erg sec“^ cm“^ after stacking twelve sources) 


is more realistic in terms of current surveys. Even the brightest 
galaxy (LI4) would hardly be detectable individually. In a sim¬ 
ulated z - 5.7 survey with depth of E(Ly(r) = 3E-18 erg sec“^ 
cm^ (~0.9E-18 erg sec“^ cm“^ after stacking twelve sources), 
the Lya surface brightness profile reached the level of the back¬ 
ground noise at a radius of 5 kpc. 

In Eig. we compare the rest-frame UV and Lya line sur¬ 
face brightness. The figure shows that the UV profiles for the 
LARS-LAEs and LARSaverage stacks are very similar. They are 
also steeper than any other Lya-stack at both redshifts, show¬ 
ing Sersic indices ^ ~ 9 and ^~llatz~2 and z - 5.7, 
respectively). The consequence of this is that the sizes of the 
Lya-stacks are generally larger than the continuum ones (Paper 
0, Table [^. The profile of the LARSaverage stack (composed 
of both EW(Lya) > 20 A and EW(Lya) < 20 A sources) was 
also shallower and more extended than that of the LARS-LAEs; 
i.e. the Lya emission of the LARS-LAE stack was more com¬ 
pact. As we described in Section |3.2[ there are no unique con¬ 
ditions for Lya photons to escape or unique morphologies. The 
conditions that make a galaxy a Lya emitter (mainly Lya flux 
concentrated around rest-frame UV bright star-forming regions) 
produce a consistent surface brightness profile, while variations 
in the dust and HI contents and distributions, which made Lya 
photons eventually escape along the line of sight, produce a more 
complex, patchy, and extended emission. 

5.3. Ground-based PSF convolution and stacking of LARS at 
high redshift 

To be able to properly compare with current observations at 
z ~ 2 and z - 5.7, we convolved LARS images with 
ground-based PSE (Eig. [n].. As the ground-based PSE is 
much larger than the HST one, we applied a Gaussian ker¬ 
nel with cr = PSE(pixel)/2.3548, using the GAUSS task in 
iraf .image s.im filter. [Eeldmeier et al.| ( |2013| ) data at z - 2.07 
exhibited a 1.4" PSE in narrow-band observation and 1" in the 
broad band. SI 1 observations at z - 2.64 were performed in sim¬ 
ilar weather conditions (PSE=0.8"-1.2"). The measurements re¬ 
ported in M14 were performed on frames smoothed to a EWHM 
= 1.32". To be able to compare with their observations at both 
z - 2.2 and z - 5.7, we convolved line and continuum images 
with that same Gaussian kernel, before stacking. This convolu¬ 
tion produced a profile close to Gaussian in the stacked images 
of both continuum and Lya; the central peak was suppressed 
and, as a consequence, the profile became shallower. However, 
the stacked Lya profile remained more extended than that of the 
rest-frame UV (Table [T^. 

In the studies mentioned earlier the surface brightness of the 
circumgalactic medium (CGM) was estimated by fitting an ex¬ 
ponential curve to radial profiles that excluded the first 2" (~ 17 
kpc at z ~ 2 and ~ 11 kpc at z ~ 5.7). Inside the first 2", emission 
from the interstellar medium (ISM) dominates over that from the 
CGM. Moreover, the ground-based PSE will most likely domi¬ 
nate on these (and smaller) scales. 

A fair comparison with high-z observations should have been 
performed to radii larger than those allowed by our HST obser¬ 
vations, in which we mainly detected Lya photons coming also 
from the interstellar media. Looking at Eig. 9 of SI 1 and Eig. 3 
of Ml4, we derived the high-z Lya profiles on scales of 10-20 
(z ~ 2) and 5-15 (z - 5.7) kpc, resembling the profiles typi¬ 
cal of high-z interstellar media, with some contribution from the 
PSE. We show them in Eig. labelled as “Steidel ISM” and 
“Momose ISM”. The surface brightness of M14 stacked sam- 
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L. Guaita: LARS at low and high z 
Table 13: Size of stacked Lyof and continuum images 


LARS stack 

Ex 

^circ 

rP20‘" 

rp20«'-^ 

j.p20fninA 

r20 

Lo 

TSO 


kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

z ~ 2 








F(Lya) = 5E - 19 








LARS-LAE 

5.26 

6.67 

6.00 

3.34 

1.33 

2.67 

5.34 

LARSaverage 

8.05 

> 12 

> 12 

8.01 

2.84 

7.84 

> 12 

F(LycK) = 3E - 18 








LARS-LAE 

3.55 

6.84 

5.00 

3.67 

1.50 

2.67 

5.67 

LARSaverage 

4.84 

> 12 

> 12 

7.00 

2.84 

7.34 

> 12 

E(LycK) = 2E - 11 








LARS-LAE 

2.16 

4.00 

5.34 

3.67 

1.17 

2.50 

4.00 

LARSaverage 

2.54 

7.00 

9.01 

4.00 

1.83 

3.84 

6.84 

z^5.7 








F{Lya) = 5E - 19 








LARS-LAE 

2.40 

4.70 

6.34 

3.52 

1.17 

2.35 

4.35 

LARSaverage 

2.72 

11.98 

11.75 

5.87 

2.11 

5.17 

10.34 

E(Lycr) = 3E - 18 








LARS-LAE 

0.81 

4.35 

3.52 

2.82 

1.29 

2.47 

4.11 

LARSaverage 

0.82 

5.05 

3.29 

4.70 

1.76 

3.17 

5.17 

z ~ 2 








^rest-UV = 30 








LARS-LAE 

2.28 

2.17 

1.33 

1.33 

0.50 

1.00 

1.83 

LARSaverage 

3.09 

4.34 

3.00 

1.67 

0.67 

1.50 

3.50 

z^5.7 

^rest-UV = 30 








LARS-LAE 

1.28 

1.88 

1.41 

0.94 

0.35 

0.82 

1.53 

LARSaverage 

1.45 

3.76 

1.64 

1.17 

0.59 

1.29 

3.29 


Notes. Size measurements from elliptical and circular aperture photometry. For three(two) simulated survey depths, we reported the Lycr emission 
sizes of the two stacks (LARS-LAE, LARSaverage) at z ~ 2(z - 5.7). We also reported the continuum emission sizes of the two stacks in the 
lower part of the table at a lOcr detection limit of mrest-uv - 30.0. F(Lycr) is given in erg sec"^ cm“^. 
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Fig. 15: Lycr surface brightness profiles of the two stacks: LARS-LAE {upper) and LARSaverage {lower, see text for details). Left and right 
panels show the profiles at z ~ 2 and 5.7, respectively, corresponding to three and two lOcr detection limits (blue circles, stars, and triangles). 
The cyan solid line represents the LARS-LAE profile of the deepest simulated survey and it is shown in the lower panels for comparison. The 
difference in the surface brightnesses in the left and right panels is the result of surface brightness dimming. The error bars correspond to the 
shallowest survey in each panel; for the other depths they are usually smaller than the symbols. The large error bars at sma > 5 kpc are produced 
by background noise. 


pie at z - 2.2 was more than one order of magnitude fainter 
than Steidel’s and ours. LARS surface brightness profile was 
as bright as Steidel’s one and steeper than Sll and M14 sam¬ 
ples. The characteristic scale of a fitted exponential profile, I(r) 
cx exp(-r/r„) was r„ = 4, almost half the value we calculated 


for Sll and M14 ISM profiles. At z ~ 5.7 we obtained a similar 
slope profile as in Ml4, with r„ ~ 2.7. 
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L. Guaita: LARS at low and high z 


F(LyQ)=5E-19, z^2 




0 LARS-LAE(LyQ) 

A LARSaverage(Lyo) 
'S' LARS-LAE(UV) 
i LARSaverage(UV) 

t t 


10-24 r* ■ ■ ' -^^^- 

0 2 4 6 8 10 12 




Fig. 16: Left’. Surface brightness profiles for Lycr (units of erg sec"^ cm“^ kpc"^) at z ~ 2 {upper panel). LARS-LAEs are shown with the cyan 
circles and line, and LARSaverage with (blue triangles). Rest UV profiles (units of erg sec“^ cm“^ kpc“^ A“^) are shown for the LARS-LAE 
sample (green stars) and LARSaverage (green line with errorbars). The lower panel shows the equivalent profiles at z = 5.7. Right. Erom the 
upper left to the lower right. Lycr for the LARS-LAE and LARSaverage stacks; rest-frame UV continuum for LARS-LAEs and LARSaverage. 
The logarithmic colour scale is adjusted to match the background level. The elongated structures in the LARSaverage continuum frame come from 
L09 (vertical) and Lll (diagonal). The whitish structures in the LARSaverage Lycr come from the same galaxies, which are weak Lycr emitters 
but show strong absorption centrally. The cut-out size is about 15 kpc. 




Eig. 17: Left: Lycr surface brightness, I(r), of the z ~ 2 LARS-LAE and LARSaverage stacks, before convolution with ground-based PSE (cyan 
circles and blue triangles). The PSE-convolved LARS-LAE stack is shown with the cyan line. The ground-based PSE profile from Eeldmeier| 
et al.|(|2013|) is shown as a black line. The blue and red dashed-dotted lines represent the circumgalactic medium, as fit by |Steidel et al.| ( ^ll| ) and 
Momose e t al.|(|201 4|) at z - 2.2, we e xtrapolated inwards. T he blue and red solid lines represent the I(r) fit on scales of 10-20 kpc derived from Eig. 


9 of |Steidel et al.|j201l| and Eig. 3 of 


Momose et al. 


__ Middle: Luminosity per kpc^ for our PSE-convolved LARS-LAE stack in Lycr (cyan 

line ), for the ground-based PSE profile from |Feldmeier et al.| ( |2013| ) (black line), for the 10-20-kpc-scale fit from |Steidel et al.|(2011| ) (blue line) 
and [Momose et al.| < |2014| ) (red line). Right: Lycr surface brightness of the z - 5.7 LARS-LAE and LARSaverage stacks before convolving with 
ground-based PSF (cyan circles and blue triangle s) together with that o f the PSE-convolved LARS-LAE stack (cyan line). The red dashed-dotted 
line represents the circumgalac tic medium fit by [Momose et al.| ( |2014| ) at z - 5.7, we extended inwards. The red line shows the fit on scales of 
5-15 kpc derived from Eig. 3 of Momose et al. ( [2014[ ). 


6. Discussion 

In Sec. we have presented the morphological properties of 
LARS galaxies; in Sec.[^ we have described the results of the 
test of simulating LARS galaxies at z ~ 2, following the meth¬ 
ods explained in Sec[^ in Sec[^ we have performed the stacking 


of individual LARS galaxy frames to simulate their typical ex¬ 
tended Lya emission at z ~ 2 and z - 5.7. Here we discuss the 
main results. 
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L. Guaita: LARS at low and high z 

Table 14: Size of the LARS-LAE stack convolved with ground-based PSF 


LARS stack 

^circ 

rP2(y" 

rp2o^'^ 


r20 

Lo 

TSO 


kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

z ~ 2 








F(Lya) = 5E-19 








LARS-LAE 

5.26 

6.67 

6.00 

3.34 

1.33 

2.67 

5.34 

Eeldmeier2013, r„ 






5.4 


LARS-LAE,PSE,E 

9.40 

9.51 

9.34 

6.00 

3.00 

5.00 

7.84 

Momose2014, r„ 






7.9 


LARS-LAE,PSE,M 

9.21 

9.34 

9.01 

5.67 

2.84 

4.84 

7.67 

Steidel2011,r“^“'’’ 






28.4 


^rest-UV = 30 








LARS-LAE 

2.28 

2.17 

1.33 

1.33 

0.50 

1.00 

1.83 

LARS-LAE,PSE,E 

5.69 

5.50 

5.34 

> 12 

2.00 

3.34 

4.67 

LARS-LAE,PSE,M 

6.40 

6.67 

6.67 

4.34 

2.50 

4.00 

5.67 

Steidel2011,r^"^^^”^^ 






2.9 


z^5.7 








E(LyQf) = 5E-19 








LARS-LAE 

2.40 

4.70 

6.34 

3.52 

1.17 

2.35 

4.35 

Momose2014, r„ 






5.9 


LARS-LAE,PSE,M 

6.60 

9.28 

7.99 

4.70 

2.35 

4.11 

7.28 

^rest-UV = 30 








LARS-LAE 

1.28 

1.88 

1.41 

0.94 

0.35 

0.82 

1.53 

LARS-LAE,PSE,M 

5.59 

5.99 

7.52 

4.46 

2.00 

3.29 

4.93 


Notes. Sizes for the LARS-LAE stack before and after convolving with |Eeldmeier et al^ ( |2013| ) (PSE,E) and |Momose et al.|f2014| ) (PSF,M) ground- 
based PSF. The values correspond to the Lycr and rest-frame UV continuum images. F or reference we report the e-foldi ng parameters of the best-fit 
exponential profile, I(r), of the observed surface brightness by |Feldmeier et al.| ( |2013| ), by |Momose et al.| ( [2014| ), and by |Steidel et al.| ( [20TT) , where 
I(r) oc exp(-r/r„) and r > 10 kpc. 


6.1. Local Universe LARS 


As described in Paper II, LARS galaxies are irregular, star¬ 
forming galaxies. Compared to the non-LAE LARS galaxies, the 
LARS-LAEs are younger and are characterized by lower star- 
formation rates, by Lycr escape fractions larger than 10% (except 
LOS), by low dust content in terms of the ratio Ha/Hfi, and by 
lower masses. Eigs|^ andtell us that the LARS galaxies have 
sizes, stellar masses, and rest-frame absolute magnitudes similar 
to those of Lyman break analogues, and that they are also com¬ 
parable to 2 < z < 3 star-forming galaxies. The stellar masses 
of the LARS galaxies tend to be larger than those estimated for 
LAEs at z ~ 2 - 3. Even if most of the LARS-LAEs have Lycr lu¬ 
minosities twice fainter than those of high-z hya emitters, they 
have stellar masses and sizes compara ble to those of the sub- 
sample of the most massive LAEs from Guaita et ( |20 11 1 ) and 
[Bond et ar] ( |2012| ). Therefore, LARS galaxies are analogues of 
10^-10^^ Mo high-z star-forming galaxies; they also share some 
properties with the most massive Lya emitters at high z. 

The ratio between half-light radii estimated in the rest-frame 
optical and rest-frame UV have been adopted in the literature 
to identify the presence of multiple stellar populations. Within 
LARS, this ratio varies from 0.8 to 3, with a median value larger 
than 1 (~ 1.4). This shows that the young stellar populations are 
more localized than the old ones, consistent with the findings for 
z ~ 2 SMGs ( [Swinbank et al.|2010| ), sBz K (|Yuma et al.|2012| ), 
and star-forming galaxies at 1.4 < z < 2.9 ( [Bond et al.|2011| ). 

The optical morphologies of LARS are typical of irregular or 
starburst galaxies and merging systems. They vary from the typ¬ 
ical of early-stage mergers to closely-gathered clumps of intense 
star-formation. G and M20 are smaller when measured in hya 
than in the rest-frame UV. However, they are comparable when 
measured in the rest-frame optical and in Uya (Eig.[^. This in¬ 
dicates quantitatively that the Uya emission of LARS galaxies 
is generally characterized by one component in a structure that 
tends to be more extended in Uya than in the bright UV contin¬ 
uum and seems to follow the entire galaxy stellar populations. 


Also, the LARS-LAEs are the ones characterized by the high¬ 
est concentration, lowest asymmetry and lowest dumpiness in 
\jya. A \jya emitting galaxy with these properties would eas¬ 
ily satisfy the EWCLya) requirement adopted at high redshift to 
identify Uya emitters. 


With the aim of identifying any continuum morphological 
property that characterizes Uya emitters and star-forming galax¬ 
ies, we study the correlations between morphological and physi¬ 
cal parameters in Table The LARS galaxies characterized by 
larger EW(LyQf) tend to be smaller in the rest-frame UV and op¬ 
tical than the other galaxies. They also present a more symmetric 
Lya emission. The galaxies which are younger (Age <10 Myr), 
less massive (M=, < 10^^ Mo), and characterized by larger spe¬ 
cific star-formation rate (dust-corrected sSERuv > -9.3 yr“^) 
tend to present larger Gini coefficient for a fixed M20 value. 
These are mainly the LARS-LAEs where the merging compo¬ 


nents appear separated, like in early-stage mergers. In Eig. 18 


we show the G vs M20 diagram, in which the colour scale cor¬ 
responds to the integrated physical parameters (Paper II). 


In Lya, LARS-LAEs tend to be composed of one bulge¬ 
like component and to harbour lower S and ellipticity. The other 
LARS galaxies are composed of patchy Lya emission (see also 
Appendix B). Some LARS galaxies with large E(B-V)neb tend 
to be characterized by the lowest M20 in the rest-frame opti¬ 
cal, typical of a bulge-like morphology. According to these re¬ 
sults, it seems that early-stage mergers could be characterized 
by younger stellar populations and symmetrical, somewhat ho¬ 
mogeneous Lya emissions, which could satisfy high-z LAE se¬ 
lections. Late-stage mergers could instead be characterized by 
turbulent star-formation episodes and, as a consequence, patchy 
Lya emissions, just like the ones observed within one galaxy 
with numerous star-forming regions. This may indicate that the 
clumpy Lyg emitting galaxies observed at high reds hift ( [Bond 
et al.|2Q12}|Gronwall et al.|2QTH|Shibuya et aL|2014|) could also 
be experiencing early-stage merging events ( jCooke et al.|20I0| ). 

Of the twelve LARS galaxies, whose continuum morphologies 
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L. Guaita: LARS at low and high z 
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Fig. 18: vs M20 measured in the rest-frame optical (left panels) and Lya (right panel) of the original LARS galaxies (squares). LARS- 

LAEs are indicated with open circles. In each raw the colour scale corresponds to an integrated physical property derived in Paper II (vertical 
colour bar). From the top raw the integrated physical properties are stellar mass, Lycr escape fraction, dust-corrected specific star formation rate, 
nebular reddening, and age. The dashed lines indicate the regions of separation between merging system, normal galaxies, and bulge-dominated 
systems as presented in Figs.[T3|and[T4 


are consistent with them being merging systems, six are LAEs. 
This does not mean that every observed LAE must necessarily 
be a merging system ( [Shibuya et al.|[MT^ |Law et al. |2012 ). 
However, we can not verify the claim by [Shibuya et al. ^20T^ 
that mergers are rare in LAEs with EW(Lya) > 100 A, because 
LARS galaxies all present lower Lya equivalent widths. 


6.2. LARS galaxies simulated to be atz- 2 

To investigate the detectability of LARS galaxies and haloes at 
high redshift, we performed simulations. The LARS galaxies 
were resampled at z ~ 2 as described in Sec 2.3. We find that 
the effects of pixel resampling and simulated survey noise are 
dependent on the irregular structures of each individual galaxy. 


Eirst of all we defined the detection limits in the continuum 
and line images (Table |^. In a survey shallower than the detec¬ 
tion limit, the background noise dominates the continuum and 
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L. Guaita: LARS at low and high z 

Table 15: Spearman coefficient and probability of the correlations between morphological and physical parameters 


parameter 

rP20*" 

QSB-rPm 

M20 

QminA 

A 

S 

ell 

^(Lya/rest-UV) 

rest-UV 









EW(LyQr) 

(-0.64,0.03) 

(0.25,0.44) 

(0.49,0.11) 

(-0.20,0.53) 

(0.22,0.50) 

(-0.29,0.37) 

(-0.11,0.74) 


fesc(Lycr) 

(-0.22,0.48) 

(0.24,0.45) 

(0.43,0.16) 

(0.09,0.77) 

(0.38,0.22) 

(0.01,0.98) 

(0.18,0.59) 


Age 

(0.52,0.08) 

(0.01,0.98) 

(-0.41,0.18) 

(0.11,0.75) 

(0.20,0.54) 

(0.68,0.01) 

(0.30,0.34) 


sSFRuv_corr 

(-0.57,0.05) 

(0.07,0.84) 

(0.25,0.43) 

(0.00,0.99) 

(-0.30,0.34) 

(-0.74,0.01) 

(-0.16,0.61) 


M, 

(0.74,0.01) 

(0.06,0.86) 

(-0.33,0.30) 

(-0.02,0.96) 

(0.09,0.78) 

(0.63,0.03) 

(0.16,0.62) 


E(B-V)neb 

(0.01,0.97) 

(-0.44,0.15) 

(-0.20,0.53) 

(-0.23,0.48) 

(-0.31,0.32) 

(-0.13,0.69) 

(-0.42,0.17) 


rest-optical 









EW(LyQr) 

(-0.60,0.04) 

(0.24,0.46) 

(-0.03,0.93) 

(-0.57,0.05) 

(-0.13,0.70) 

(-0.36,0.25) 

(-0.18,0.57) 


fesc(Lycr) 

(-0.38,0.22) 

(0.11,0.73) 

(0.24,0.44) 

(-0.31,0.33) 

(0.01,0.98) 

(0.03,0.93) 

(0.24,0.46) 


Age 

(0.58,0.05) 

(-0.62,0.03) 

(0.05,0.88) 

(-0.17,0.60) 

(-0.39,0.21) 

(0.26,0.42) 

(0.01,0.98) 


sSERuv_corr 

(-0.69,0.01) 

(0.62,0.03) 

(-0.16,0.62) 

(0.08,0.81) 

(0.50,0.10) 

(-0.27,0.40) 

(-0.01,0.98) 


M, 

(0.66,0.02) 

(-0.55,0.07) 

(0.24,0.46) 

(0.21,0.51) 

(-0.13,0.70) 

(0.36,0.25) 

(0.14,0.66) 


E(B-V)neb 

(0.08,0.81) 

(-0.32,0.31) 

(-0.42,0.17) 

(0.20,0.54) 

(-0.10,0.75) 

(-0.18,0.58) 

(-0.41,0.18) 


Lycr 









EW(LyQr) 

(-0.22,0.53) 

(0.51,0.09) 

(0.03,0.93) 

(0.59,0.04) 

(-0.70,0.01) 

(-0.61,0.04) 

(-0.42,0.17) 

(0.25,0.49) 

fesc(Lycr) 

(-0.52,0.13) 

(0.20,0.53) 

(0.42,0.17) 

(0.18,0.57) 

(-0.67,0.02) 

(-0.76,0.00) 

(-0.19,0.56) 

(-0.22,0.53) 

Age 

(-0.01,0.99) 

(0.60,0.04) 

(0.22,0.50) 

(-0.32,0.31) 

(0.61,0.04) 

(0.84,0.00) 

(0.57,0.05) 

(-0.24,0.51) 

sSERuv_corr 

(0.02,0.96) 

(0.65,0.02) 

(-0.38,0.23) 

(0.27,0.40) 

(-0.56,0.06) 

(-0.90,0.00) 

(-0.54,0.07) 

(0.49,0.15) 

M, 

(-0.03,0.93) 

(-0.45,0.14) 

(0.31,0.33) 

(-0.50,0.10) 

(0.76,0.00) 

(0.80,0.00) 

(0.56,0.06) 

(-0.47,0.17) 

E(B-V)neb 

(0.53,0.12) 

(-0.19,0.55) 

(-0.34,0.29) 

(0.02,0.95) 

(0.52,0.08) 

(0.57,0.05) 

(-0.06,0.85) 

(0.31,0.38) 


Notes. The number pairs in the table correspond to the Spearman coefficient, r, and probability, p, for the correlations between physical (column, y) 
and morphological (row, x) parameters, measured in the rest-frame UV, optical, and Lycr images. For instance, for and y=sSFRuv_corr, 

the Spearman(x,y) calculation tells us that the probability that x and y are uncorrelated datasets is only p=3% and there is a positive (r = 0.62) 
correlation. Bold face indicates the most significant correlations. 


line surface brightness profiles. Interestingly all LARS galax¬ 
ies would be detected in the rest-frame UV in surveys like 
HUDF09 and GOODS, and at the rest-frame optical in a survey 
like CANDELS (deep), if located at z ~ 2. However, even more 
interesting is that Lya emission extended up to 5 kpc would be 
visible in 70% of the sample at a lOcr depth of 3E-18 erg sec“^ 
cm“^ (Table and Appendix C). The LARS galaxies charac¬ 
terized by the faintest integrated Lya flux show an even more 
extended emission. This should be taken into account when de¬ 
signing Lya spectroscopic observations (see also the discussion 
in Paper I), which aim to detect as much of the Lya flux as 
possible. In the current Lya surveys (typically characterized by 
E(LyQf) > 2E-17 erg sec“^ cm“^) only L14, the brightest LARS- 
LAE, would be detected. Even if the Lya size can be as large as 
5 times the UV, the extremely detailed Lya structures (Eigs. 
and [T]) could only be identified in the original HST-resolution 
images. 


Second of all we quantified the morphological parameters of 
high-z LARS galaxies. The z2LARS-LAEs tend to be smaller 
in the rest-frame UV, optical, and Her than the other LARS 
galaxies. This is consistent with the observed high-z findings 
dGronwall et aLllMTT] [Bond et al.|[MT2l ). Also, the z2LARS- 
LAEs present a large range of asymmetry values in UV and op¬ 
tical, but do show symmetric morphologies of Lya emission. 

The low surface brightness structures tend to disappear 
within the background in shallow surveys, making a clumpier 
and sometimes more elliptical galaxy in the continuum. As a 
consequence the ellipticity values we measured in the rest-frame 
UV and optical of z2LARS-LAEs increase as the survey depth 
decreased. Also, the continuum dumpiness tends to significantly 
increase. Some of the high-z observations of large elli pticity and 
clumpy systems could be expla ined in terms of depth ( |GronwalI| 
|et al.||20Il Bond et al. |20I2| ), but we cannot explain the de¬ 
crease in the merg er fraction and elliptic ity for EW(LyQf) >100 
A LAEs found by Shibuya et al. ( 2014| . However, the dumpi¬ 
ness measured i n Ly a als o increases for decreasing survey depth 
(Eigs. |B.l] |B.2| and |B.3| ), making the integrated EW more dif¬ 


ficult to estimate. This happens for the LARS-LAEs as well, 
which eventually become close to impossible to detect as Lya 
emitters. In our sample LOS is the LARS galaxy with the low¬ 
est rest-frame optical ellipticity in the original, high-z simu¬ 
lated, and shallow-survey images. It is a massive face-on irregu¬ 
lar galaxy in the rest-frame optical, wit h mul tiple star-formation 
clumps seen in the rest-frame UV (Eig. |B.2| ), but it is character¬ 
ized by the lowest EW(LyQf) among the LARS-LAEs. 


Continuum G and M20 values are preserved after pixel re¬ 
sampling and adding noise (see also Appendix A). Therefore, we 
adopted these two parameters for characterizing LARS galaxies 
at high redshift (Eig. [^. LARS galaxies have a morphology 
consistent with merging systems even when simulated at high z. 
Some z2LARS-LAEs have both rest-frame optical and UV mor¬ 
phologies consistent with being mergers (Eig. El*. 


The asymmetry we estimated decreases after adding noise. A 
merger information calculated just adopting a large-asymmetry 
criterion could lead to a mis-interpretation of our sample in a 
shallow survey. This w as also noticed in [Shibuya et al. ( |20I4| ); 


Gronwall et al. (2011) for low signal-to-noise sources. 


In general the Lya morphology tends to be significantly af¬ 
fected in shallow surveys, because Lya detailed structures are 
very sensitive to depth and resolution. 


6.3. Lya haloes of LARS galaxies and their implications 

There is still open debate about the conditions necessary for the 
formation of Lya haloes in high-z galaxies, but it seems that HI 
scattering is the main factor at the scale of LARS galaxies (Paper 
II). At high redshift, Lya halo studies have been performed in 
stacked data obtained from various samples of galaxies. To be 
able to compare with high-z results, we simulated how the Lya 
haloes of LARS galaxies would appear at z ~ 2 and z - 5.7. 
This was performed by stacking subsamples of our galaxies and 
assuming a range of survey depths. In this test we have the ad¬ 
vantage of knowing the original morphology and halo profile. 
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By simply examining the RGB mosaics in Fig. 1, we see 
that stacking images with very different irregular structures can 
be very complicated, but it is commonly done at high redshift 
to increase the signal-to-noise. Background noise, survey depth, 
and ground-based Point Spread Function are the primary limits 
of the detection of Lya haloes. We find that a depth compara¬ 
ble to the Ml4 survey is ideal to r ecover these haloes at both 
z ~ 2 and z - 5.7 (see Sec. IMI Fig.[Tg. The stacked Lya 
surface brightness profile of the LARS-LAEs is peaked in the 
centre, whereas the LARSaverage stack contains contributions 
from more diverse (asymmetric, patchy, and eventually more ex¬ 
tended) Lya morphologies; the result is a more irregular and 
extended profile. It seems reasonable that the conditions under 
which a galaxy may emit Lya (mainly Lya flux concentrated 
around the UV-bright knots of star formation), would produce a 
consistent surface brightness profile. However, variations in the 
dust and HI contents and distributions, which may cause Lya 
photons escape along the line of sight (EW(LyQf) > 0 A, but not 
necessarily > 20 A), would produce a more complex, patchy, 
and extended emission. This is seen in the LARSaverage stack 
and also in the stack observed in high-z LBGs. 

Ground-based observations of high-z galaxies allow us to 
construct Lya profiles which extend much farther from the cen¬ 
tre (see Sll and Ml4) than low-z galaxies observed with the 
Hubble Space Telescope. The continuum subtraction procedure 
( [Hayes et al.||2009| ), described in Paper II to isolate Lya emis¬ 
sion, takes advantage of the HST resolution to provide a detailed 
Lya mapping within and just outside the interstellar medium. 
Within the first 10 kpc Lya scattering already begins to produce 
halo-like structures in LARS galaxies, while it is at radii above 
10-20 kpc that high-z studies are performed because no spatial 
information is available inside the PSF. A fair comparison be¬ 
tween low and high-z observations could only be made by ap¬ 
plying the same procedure to extract Lya from the rest-frame 
UV emission at high redshift and by investigating the medium 
on the same scales and the same PSF conditions. 

We convolved our HST images with typical ground-based 
PSFs and compared our smoothed Lya profile with ground- 
based observations at comparable scales. Even in this case the 
z ~ 2 LARS Lya stacked profile is steeper than the one derived 
for the complete sample by M14 and for the LAE-only subsam¬ 
ple of Sll. In Feldmeier’s survey, the LARS Lya profile would 
be indistinguishable from the ground-based PSF on scales larger 
than 6 kpc. At z ~ 5.7, LARS Lya profile is as steep as the one 
from Ml4, which is brighter than the profile they obtained at 
z ~ 2. 

With LARS we cannot probe as large scale as the current 
high-z observations do, due to smaller field-of-view of HST. 
However, we could still expect some differences between LARS 
and Feldmeier’s sample because of observational depth, and be¬ 
tween LARS and the S11 sample because LARS are not located 
in overdense regions. Also, we may expect differences between 
LARS and M14 sample at large scale due to the difference in 
physical properties of the two samples. 

The z - 2.2 LAE sample studied by Momose et al. is char¬ 
acterized by dust reddening (E(B-V) < 0.1), metallicity (Z ~ 0.2 
Zo), and stellar mass (M=, < 10^^ M©) lower than LARS. Since 
Lya photons are sensitive to the presence of dust grains and to 
the scattering on neutral hydrogen (e.g.. Paper II), Lya morphol¬ 
ogy (extension and features) is expec ted to depend on the larger 
dust and HI contents (Paper III and [Rivera-Thorsen & Hayesj 
2015| ). The dust grains, able to absorb Lya photons close to the 
knots of star-formation, also prevent their escape at large scales. 


where HI scattering plays the role of making haloes. There are 
other phenomena, such as gas kinematics, we are investigating 
within the LARS survey, which could favour the escape of Lya 
photons from the HII regions and ultimately allow the forma¬ 
tion of Lya haloes. [Rivera-Thorsen & Hayesj ( [2015[ ), Duval et 
al., and Orlitova et al. (in prep) are dedicated to studying the gas 
kinematics in LARS galaxies from HST spectroscopy. 

7. Conclusions 

This paper is number IV of a series. In this work we have charac¬ 
terized and quantitatively studied the morphology of a sample of 
starburst galaxies at z < 0.2: the Lyman alpha reference sample, 
LARS. 

- LARS galaxies have continuum sizes and stellar masses sim¬ 
ilar to those of local Lyman break analogues and 2 < z < 
3 star-forming galaxies. The stellar mass and luminosities 
also match the two samples, respectively. Therefore, LARS 
galaxies can be studied as a reference of 10^-10^^ M© high-z 
star-forming galaxies; they also share some properties with 
the most massive Lya emitters at high redshift. 

- The rest-frame optical morphology of LARS galaxies is the 
typical of merging systems. This is also valid for the LAEs 
within LARS. 

- For the first time we were able to quantify the morphology of 
Lya emission. LARS-LAEs are on average characterized by 
more concentrated and symmetrical, while LARS non-LAEs 
can present patchier and irregular Lya emissions. LARS- 
LAEs are more compact in Lya, even when regridded to high 
redshift. 

- We have simulated LARS galaxies at high redshift and ex¬ 
plored their detection: all LARS galaxies would be detected 
in the continuum in current deep surveys, but they would not 
be easily detected in the current Lya surveys at z ~ 2. 

- In a shallow survey, it is the morphology of Lya that is 
most affected by background noise, because the detailed Lya 
structures strongly depend on depth and resolution. This may 
affect high-z Lya observations. 

- The measured ellipticity and dumpiness tend to increase in 
shallow surveys for most of the LARS galaxies. This could 
explain some of the high-z observations of large ellipticity 
and clumpy systems in LAE samples. 

- We stacked the Lya images of LARS galaxies simulated at 
high redshift. The LARS-LAE stack is peaked in the centre, 
whereas the LARSaverage stack contains contributions from 
more diverse Lya morphologies resulting in a more irregu¬ 
lar and extended profile. Variations in the dust and HI con¬ 
tents and distributions may produce more complex, patchy, 
and extended emission like the one seen in the LARSaverage 
stack and also in the stack observed in high-z LBGs. 

- The Lya haloes we study in LARS galaxies probe much- 
smaller-scale media than high-z observations. We find that 
LARS-halo-profile slope is steeper than z ~ 2 and as steep as 
with z - 5.7 observations at ~10-kpc scales, after applying 
ground-based PSE. 

A sample like LARS, at slightly larger redshift, could allow 
studying circumgalactic medium-scale haloes and relating them 
to other galaxy properties. Physical properties (as presented in 
Paper II), Hl-mass estimations (from Paper III), and kinemat¬ 
ics (analysed in [Rivera-Thorsen & Hayes|2Q15[ Duval et al., and 
Orlitova et al.) already helped in clarifying the Lya-photon prop¬ 
agation within the interstellar medium. With the large-scale halo 
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information we would be able to also investigate the mechanisms 
that transport Lya photons from the interstellar to the circum- 
galactic medium. 
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Table 5: Size of high-z simulated LARS galaxies 
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Table 6: Size of high-z simulated LARS galaxies 
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1.93 

0.21 

1.71 

0.33 

1.74 

0.15 

8E-18 

1.96 

0.37 

1.76 

0.47 

1.86 

0.23 

Ha 

2E-19 

0.92 

0.06 

1.00 

0.20 

1.05 

0.13 

6E-19 

0.93 

0.07 

0.93 

0.18 

1.13 

0.15 

lE-18 

0.96 

0.07 

0.98 

0.31 

1.15 

0.14 

3E-18 

1.03 

0.11 

1.03 

0.26 

1.24 

0.15 

L07 

rest-UV 







30 

1.52 

0.18 

1.00 

0.34 

0.67 

0.17 

29 

1.57 

0.20 

1.05 

0.36 

0.67 

0.17 

28 

1.57 

0.24 

1.12 

0.39 

0.72 

0.18 

rest-optical 

29 

4.14 

0.25 

3.34 

0.34 

1.72 

0.19 

28 

4.02 

0.21 

3.39 

0.38 

1.67 

0.17 

27 

4.21 

0.36 

3.31 

0.78 

1.73 

0.20 

26 

3.99 

0.65 

3.00 

1.12 

1.62 

0.27 

Lycr 

5E-19 

4.54 

1.30 

3.97 

0.39 

4.15 

2.19 

3E-18 

3.59 

0.59 

3.85 

0.84 

3.22 

0.39 

8E-18 

3.12 

0.72 

3.45 

1.28 

3.08 

0.58 

Ha 

6E-19 

0.32 

0.04 

0.30 

0.10 

0.48 

0.10 

wE-19 

0.32 

0.04 

0.30 

0.10 

0.48 

0.10 

lE-18 

0.32 

0.04 

0.31 

0.10 

0.47 

0.10 

3E-18 

0.34 

0.08 

0.33 

0.12 

0.48 

0.10 

L08 

rest-UV 







30 

7.94 

2.62 

6.55 

1.13 

4.40 

1.69 

29 

7.57 

2.13 

5.92 

2.35 

3.45 

1.08 

28 

6.67 

2.05 

5.85 

2.75 

3.71 

0.99 

rest-optical 

29 

10.32 

1.55 

8.67 

0.38 

4.48 

1.40 

28 

9.86 

1.89 

8.53 

0.54 

4.33 

1.64 

27 

9.51 

1.97 

8.29 

0.69 

4.06 

1.40 

26 

9.17 

1.82 

7.67 

0.96 

3.94 

1.28 

25 

8.17 

1.28 

7.67 

0.96 

3.77 

0.27 

Lycr 

5E-19 

2.23 

0.16 

1.31 

0.37 

1.58 

0.12 

3E-18 

1.35 

0.44 

0.73 

0.28 

1.34 

0.39 

8E-18 

0.76 

1.37 

0.53 

1.59 

0.82 

1.48 

Ha 2E-19 

0.67 

0.03 

0.64 

0.08 

0.86 

0.07 

6E-19 

0.63 

0.03 

0.57 

0.05 

0.80 

0.04 

lE-18 

0.63 

0.04 

0.56 

0.05 

0.79 

0.04 

3E-18 

0.95 

0.27 

0.60 

0.11 

1.13 

0.37 


Notes, continuation of Tab. Sizes in Lya and Ha images are given as ^(Lya) = sizelre^t-uv) and ^ (Her) = size{rest^opticai) ’ where the sizes are rP20^^^ 
rP20'”'”^, and r50. ^(Lycr) is the most meaningful due to the resonant nature of Lycr. 
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L. Guaita: LARS at low and high z 
Table 7: Size of high-z simulated LARS galaxies 


continuum 

iP20“‘ 

crr?20“‘ 

rP20"''^ 

(TrP20"^'^ 

Uo 

O' r50 


kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

line 

^(rP2(y") 

cr^(rP20'’") 


cr^(rP20'”'”'^) 

^(Uo) 


L09 

rest-UV 







30 

20.34 

3.32 

2.18 

1.31 

10.34 

1.50 

29 

5.60 

2.65 

5.67 

2.57 

3.19 

1.84 

28 

5.09 

3.56 

1.41 

1.03 

2.34 

2.19 

rest-optical 

29 

21.71 

0.41 

18.01 

0.33 

12.17 

0.17 

28 

21.81 

0.80 

18.08 

0.49 

12.44 

0.23 

27 

15.72 

8.35 

10.87 

5.30 

9.00 

4.68 

Lya 

5E-19 

1.27 

0.04 

8.24 

0.77 

1.61 

0.04 

Ua 

2E-19 

0.66 

0.19 

0.19 

0.06 

0.50 

0.14 

6E-19 

0.80 

0.12 

0.17 

0.06 

0.66 

0.09 

lE-18 

0.73 

0.25 

0.19 

0.08 

0.67 

0.31 

3E-18 

0.43 

0.19 

0.13 

0.05 

0.37 

0.22 

lE-17 

0.19 

0.08 

0.10 

0.04 

0.14 

0.06 

LIO 

rest-UV 







30 

4.68 

0.43 

3.23 

0.41 

2.04 

0.18 

29 

4.64 

0.66 

3.34 

0.80 

2.10 

0.23 

28 

4.39 

1.19 

3.45 

1.56 

2.10 

0.37 

rest-optical 

29 

9.65 

0.36 

6.37 

0.37 

4.08 

0.20 

28 

9.22 

0.49 

6.16 

0.50 

3.85 

0.20 

27 

9.09 

0.75 

5.72 

0.86 

3.85 

0.27 

26 

8.28 

1.23 

5.34 

1.59 

3.76 

0.45 

hya 

5E-19 

2.16 

1.13 

2.59 

0.47 

3.08 

1.55 

Ua 

2E-19 

0.52 

0.07 

0.48 

0.09 

0.54 

0.05 

6E-19 

0.52 

0.07 

0.48 

0.09 

0.54 

0.05 

lE-18 

0.50 

0.08 

0.43 

0.12 

0.53 

0.06 

Lll 

rest-UV 







30 

28.66 

1.20 

13.18 

0.64 

15.17 

0.31 

29 

25.51 

2.83 

9.64 

3.72 

13.50 

0.72 

28 

14.59 

4.83 

6.46 

2.49 

9.15 

2.73 

rest-optical 

29 

25.29 

6.73 

17.64 

3.71 

12.49 

2.77 

28 

28.06 

1.24 

17.12 

1.17 

13.93 

0.23 

27 

24.95 

1.90 

13.64 

2.29 

13.67 

0.66 

26 

18.75 

3.33 

8.51 

2.63 

11.32 

1.53 

hya 

5E-19 

0.25 

0.77 

0.81 

2.04 

0.31 

1.17 

3E-18 

0.15 

0.08 

0.24 

0.33 

0.14 

0.12 

Ha 

2E-19 

1.06 

0.11 

0.70 

0.07 

0.92 

0.05 

6E-19 

0.81 

0.08 

0.56 

0.09 

0.85 

0.10 

lE-18 

0.90 

0.73 

0.61 

0.20 

1.12 

0.80 

3E-18 

0.58 

0.10 

0.38 

0.12 

0.66 

0.07 


Notes, continuation of Tab. [5] 
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L. Guaita: LARS at low and high z 
Table 8: Size of high-z simulated LARS galaxies 


continuum 

iP20“‘ 

kpc 

crTP20“‘ 

kpc 

kpc 

crrP20"^'^ 

kpc 

Uo 

kpc 

O-TSO 

kpc 

line 

^(rP2(y") 

cr^(rP20'’") 

^(rP20'^'"'^) 

cr^(rP20'”'”'^) 

^(Uo) 


L12 

rest-UV 

30 

3.38 

0.18 

1.67 

0.33 

1.33 

0.17 

29 

3.44 

0.20 

1.67 

0.34 

1.33 

0.17 

28 

3.55 

0.25 

1.63 

0.41 

1.32 

0.18 

27 

3.12 

0.94 

1.70 

0.68 

1.10 

0.35 

rest-optical 

29 

4.79 

0.21 

1.68 

1.02 

2.05 

0.18 

28 

4.83 

0.24 

2.25 

1.47 

2.05 

0.19 

27 

5.09 

0.30 

2.82 

1.59 

2.19 

0.22 

26 

5.66 

1.53 

2.33 

1.58 

2.66 

1.24 

Lycr 

5E-19 

4.14 

0.55 

6.87 

0.31 

5.65 

0.87 

3E-18 

3.81 

1.70 

6.55 

1.07 

5.67 

2.43 

Ua 

2E-19 

0.74 

0.17 

2.04 

0.38 

0.86 

0.41 

6E-19 

0.71 

0.04 

1.99 

0.20 

0.81 

0.08 

lE-18 

0.73 

0.04 

1.99 

0.20 

0.82 

0.09 

3E-18 

0.73 

0.07 

1.84 

0.34 

0.84 

0.09 

L13 

rest-UV 

30 

5.29 

2.55 

1.66 

0.34 

2.40 

2.26 

29 

4.56 

2.08 

1.59 

0.37 

1.80 

1.92 

28 

2.27 

0.72 

1.54 

0.38 

0.93 

0.22 

rest-optical 

29 

12.57 

2.48 

6.47 

0.37 

5.94 

0.91 

28 

13.43 

0.33 

6.35 

0.37 

6.44 

0.20 

27 

9.31 

3.13 

6.02 

0.43 

5.70 

1.71 

Lya 

5E-19 

0.62 

0.14 

2.48 

3.87 

1.11 

0.34 

3E-18 

0.84 

0.60 

1.74 

0.86 

1.16 

1.03 

8E-18 

0.96 

0.58 

1.63 

0.75 

1.36 

0.99 

Ua 

2E-19 

0.61 

0.12 

0.31 

0.05 

0.55 

0.14 

6E-19 

0.63 

0.12 

0.31 

0.05 

0.57 

0.14 

lE-18 

0.58 

0.13 

0.31 

0.05 

0.48 

0.18 

3E-18 

0.47 

0.13 

0.28 

0.06 

0.39 

0.13 

L14 

rest-UV 

30 

1.67 

0.17 

1.33 

0.33 

0.83 

0.17 

29 

1.83 

0.18 

1.33 

0.33 

0.84 

0.17 

28 

1.99 

1.18 

1.19 

0.38 

1.06 

0.99 

rest-optical 

29 

4.53 

0.18 

1.20 

0.37 

1.50 

0.17 

28 

4.55 

0.46 

1.10 

0.37 

1.56 

0.41 

27 

4.88 

2.60 

1.48 

0.59 

3.11 

2.02 

Lya 

5E-19 

3.55 

0.33 

2.50 

0.25 

3.00 

0.22 

3E-18 

3.87 

0.77 

2.49 

0.27 

3.31 

1.17 

8E-18 

3.50 

0.53 

2.45 

0.29 

3.09 

0.25 

2E-17 

3.02 

0.51 

2.44 

0.31 

2.87 

0.27 

5E-17 

2.96 

0.61 

2.44 

0.49 

2.93 

0.38 

Ha 

2E-19 

0.51 

0.04 

1.39 

0.28 

0.79 

0.13 

6E-19 

0.54 

0.04 

1.39 

0.28 

0.86 

0.12 

lE-18 

0.56 

0.04 

1.39 

0.28 

0.88 

0.12 

3E-18 

1.17 

0.60 

1.39 

0.30 

2.34 

1.57 

lE-17 

0.45 

0.07 

1.22 

0.32 

0.70 

0.12 


Notes, continuation of Tab. [5] 
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L. Guaita: LARS at low and high z 

Table 9: Morphological parameters of high-z simulated LARS galaxies 


depth 

Qrno 


QminA 

O-CminA 

M20 

crM20 

A 

crA 

S 

o-S 

ell 

crell 

LOl 

rest-UV 

30 

0.73 

0.01 

2.39 

0.00 

-1.08 

0.03 

0.36 

0.01 

0.20 

0.01 

0.56 

0.02 

29 

0.70 

0.03 

2.41 

0.12 

-1.16 

0.34 

0.28 

0.02 

0.18 

0.05 

0.50 

0.06 

28 

0.67 

0.04 

2.46 

0.45 

-1.14 

0.41 

0.17 

0.04 

0.20 

0.08 

0.43 

0.09 

rest-optical 

29 

0.60 

0.00 

2.13 

0.00 

-0.91 

0.08 

0.26 

0.01 

0.06 

0.01 

0.58 

0.01 

28 

0.61 

0.01 

2.16 

0.31 

-0.98 

0.08 

0.18 

0.02 

0.11 

0.01 

0.59 

0.02 

27 

0.61 

0.02 

2.59 

0.57 

-1.03 

0.10 

0.10 

0.03 

0.22 

0.04 

0.58 

0.03 

26 

0.60 

0.03 

2.53 

0.67 

-1.17 

0.19 

0.04 

0.02 

0.47 

0.12 

0.54 

0.06 

Lyof 

5E-19 

0.62 

0.02 

2.88 

0.18 

-1.66 

0.20 

0.28 

0.04 

0.23 

0.06 

0.30 

0.03 

3E-18 

0.57 

0.01 

2.54 

0.40 

-1.60 

0.15 

0.15 

0.02 

0.49 

0.09 

0.33 

0.09 

8E-18 

0.49 

0.02 

2.11 

0.58 

-1.61 

0.28 

0.09 

0.04 

0.70 

0.15 

0.20 

0.08 

Ua 

2E-19 

0.67 

0.01 

2.39 

0.00 

-1.10 

0.26 

0.36 

0.00 

0.20 

0.07 

0.43 

0.02 

6E-19 

0.68 

0.01 

2.39 

0.00 

-0.87 

0.16 

0.35 

0.01 

0.30 

0.06 

0.46 

0.02 

IE-18 

0.67 

0.01 

2.39 

0.00 

-0.94 

0.21 

0.34 

0.01 

0.26 

0.08 

0.44 

0.04 

3E-18 

0.64 

0.02 

2.40 

0.28 

-0.97 

0.21 

0.26 

0.02 

0.26 

0.07 

0.42 

0.05 

L02 

rest-UV 

30 

0.64 

0.02 

2.40 

0.64 

-1.41 

0.22 

0.21 

0.03 

0.36 

0.06 

0.69 

0.03 

29 

0.62 

0.03 

2.29 

0.65 

-1.29 

0.27 

0.10 

0.04 

0.39 

0.06 

0.71 

0.04 

rest-optical 

29 

0.62 

0.01 

3.68 

0.19 

-1.81 

0.08 

0.02 

0.01 

0.36 

0.04 

0.63 

0.02 

28 

0.60 

0.01 

3.38 

0.32 

-1.78 

0.12 

0.01 

0.01 

0.54 

0.08 

0.66 

0.03 

27 

0.55 

0.03 

2.67 

0.59 

-1.59 

0.19 

0.01 

0.02 

0.60 

0.13 

0.68 

0.05 

Lyof 

5E-19 

0.61 

0.01 

3.06 

0.35 

-1.25 

0.07 

0.27 

0.01 

0.35 

0.04 

0.36 

0.03 

3E-18 

0.51 

0.01 

2.84 

0.43 

-1.22 

0.21 

0.09 

0.03 

0.64 

0.07 

0.36 

0.12 

8E-18 

0.45 

0.02 

2.05 

0.52 

-1.20 

0.25 

0.05 

0.04 

0.68 

0.19 

0.39 

0.12 

Ha 

2E-19 

0.66 

0.02 

2.51 

0.34 

-1.29 

0.12 

0.34 

0.02 

0.15 

0.07 

0.28 

0.14 

6E-19 

0.60 

0.02 

3.50 

0.00 

-0.93 

0.12 

0.26 

0.03 

0.25 

0.07 

0.61 

0.04 

lE-18 

0.58 

0.02 

3.10 

0.79 

-0.96 

0.19 

0.21 

0.04 

0.26 

0.07 

0.62 

0.05 

L03 

rest-UV 

30 

0.67 

0.03 

2.54 

0.27 

-1.42 

0.29 

0.08 

0.03 

0.19 

0.12 

0.14 

0.08 

29 

0.64 

0.15 

2.82 

1.66 

-1.58 

1.49 

0.06 

0.14 

0.38 

1.02 

0.20 

0.50 

rest-optical 

29 

0.63 

0.06 

3.38 

0.10 

-2.01 

0.32 

0.31 

0.03 

0.09 

0.07 

0.46 

0.07 

28 

0.61 

0.02 

3.07 

0.28 

-2.00 

0.19 

0.10 

0.03 

0.13 

0.03 

0.46 

0.04 

27 

0.63 

0.03 

2.98 

0.36 

-1.96 

0.33 

0.06 

0.01 

0.29 

0.07 

0.34 

0.07 

26 

0.61 

0.03 

3.06 

0.34 

-2.04 

0.13 

0.03 

0.01 

0.53 

0.11 

0.27 

0.07 

25 

0.57 

0.02 

3.06 

0.34 

-1.87 

0.19 

0.03 

0.01 

0.64 

0.13 

0.37 

0.10 

Lyof 

5E-19 

0.47 

0.01 

2.70 

0.12 

-1.78 

0.15 

0.10 

0.02 

0.62 

0.03 

0.30 

0.06 

3E-18 

0.41 

0.09 

2.04 

2.74 

-0.93 

0.90 

0.02 

0.10 

0.76 

0.60 

0.50 

0.75 

Ua 

2E-19 

0.67 

0.00 

3.01 

0.00 

-1.58 

0.25 

0.25 

0.00 

0.08 

0.01 

0.21 

0.03 

6E-19 

0.66 

0.01 

3.01 

0.00 

-1.57 

0.25 

0.21 

0.01 

0.10 

0.01 

0.27 

0.03 

lE-18 

0.67 

0.01 

3.02 

0.05 

-1.56 

0.27 

0.19 

0.01 

0.11 

0.01 

0.29 

0.03 

3E-18 

0.65 

0.01 

2.99 

0.33 

-1.60 

0.26 

0.14 

0.03 

0.17 

0.04 

0.35 

0.05 


Notes. Non-parametric measurements are derived by following the equations in Appendix A. Only and are reported here, because they 
are the ones used to compare z2LARS with high-z observations. The uncertainty on the parameters is obtained from the standard deviation of 100 
realizations of a simulated survey depth. 
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L. Guaita: LARS at low and high z 

Table 10: Morphological parameters of high-z simulated LARS galaxies 


depth 

QrPlO 


QminA 

O-Cn^nA 

M20 

O-M20 

A 

crA 

S 

a-S 

ell 

crell 

L05 

rest-UV 













30 

0.70 

0.01 

3.88 

0.07 

-1.22 

0.07 

0.27 

0.01 

0.23 

0.04 

0.67 

0.01 

29 

0.70 

0.01 

3.73 

0.21 

-1.13 

0.10 

0.19 

0.02 

0.35 

0.12 

0.71 

0.02 

28 

0.70 

0.01 

3.68 

0.48 

-1.13 

0.12 

0.11 

0.03 

0.32 

0.14 

0.73 

0.03 

27 

0.68 

0.03 

2.82 

0.78 

-1.22 

0.25 

0.05 

0.04 

0.36 

0.13 

0.73 

0.05 

rest-optical 

29 

0.74 

0.01 

3.89 

0.00 

-1.42 

0.17 

0.27 

0.01 

0.29 

0.05 

0.73 

0.01 

28 

0.74 

0.01 

3.90 

0.06 

-1.40 

0.18 

0.19 

0.02 

0.27 

0.07 

0.76 

0.02 

27 

0.75 

0.02 

3.93 

0.24 

-1.28 

0.24 

0.12 

0.03 

0.25 

0.11 

0.79 

0.02 

26 

0.74 

0.03 

3.31 

1.04 

-1.26 

0.29 

0.06 

0.03 

0.33 

0.17 

0.81 

0.03 

Lyor 

5E-19 

0.58 

0.01 

1.95 

0.09 

-1.42 

0.06 

0.26 

0.00 

0.11 

0.01 

0.13 

0.03 

3E-18 

0.52 

0.02 

1.98 

0.38 

-1.48 

0.08 

0.13 

0.02 

0.39 

0.08 

0.17 

0.05 

8E-18 

0.45 

0.01 

2.18 

0.38 

-1.40 

0.13 

0.05 

0.03 

0.64 

0.11 

0.27 

0.10 

Ha 

2E-19 

0.59 

0.01 

2.21 

0.20 

-0.79 

0.01 

0.31 

0.01 

0.27 

0.01 

0.44 

0.02 

6E-19 

0.57 

0.01 

2.25 

0.19 

-0.76 

0.02 

0.26 

0.03 

0.24 

0.04 

0.52 

0.02 

lE-18 

0.56 

0.01 

2.21 

0.20 

-0.76 

0.02 

0.21 

0.04 

0.22 

0.04 

0.56 

0.03 

3E-18 

0.52 

0.02 

2.24 

0.24 

-0.75 

0.09 

0.11 

0.04 

0.27 

0.06 

0.64 

0.04 

L07 

rest-UV 













30 

0.72 

0.01 

2.39 

0.00 

-1.20 

0.11 

0.31 

0.01 

0.02 

0.05 

0.47 

0.02 

29 

0.72 

0.01 

2.39 

0.00 

-1.12 

0.23 

0.27 

0.01 

0.07 

0.10 

0.50 

0.03 

28 

0.72 

0.02 

2.43 

0.24 

-1.02 

0.29 

0.19 

0.04 

0.14 

0.12 

0.51 

0.06 

rest-optical 

29 

0.63 

0.01 

3.01 

0.03 

-1.54 

0.12 

0.23 

0.01 

0.05 

0.01 

0.59 

0.03 

28 

0.62 

0.01 

3.05 

0.13 

-1.51 

0.12 

0.15 

0.01 

0.08 

0.01 

0.59 

0.02 

27 

0.63 

0.01 

2.97 

0.38 

-1.54 

0.13 

0.07 

0.02 

0.20 

0.04 

0.61 

0.02 

26 

0.63 

0.02 

2.70 

0.65 

-1.51 

0.18 

0.04 

0.03 

0.46 

0.12 

0.55 

0.09 

hya 

5E-19 

0.64 

0.01 

3.44 

0.10 

-2.17 

0.21 

0.24 

0.01 

0.18 

0.05 

0.36 

0.04 

3E-18 

0.57 

0.01 

3.18 

0.27 

-2.18 

0.12 

0.11 

0.02 

0.49 

0.09 

0.13 

0.07 

8E-18 

0.51 

0.02 

2.85 

0.48 

-2.08 

0.16 

0.06 

0.03 

0.66 

0.12 

0.17 

0.08 

Ha 

2E-19 

0.65 

0.01 

0.88 

0.00 

-0.72 

0.06 

0.34 

0.01 

0.23 

0.05 

0.24 

0.07 

6E-19 

0.65 

0.01 

0.88 

0.00 

-0.72 

0.06 

0.34 

0.01 

0.23 

0.05 

0.24 

0.07 

lE-18 

0.66 

0.01 

0.88 

0.00 

-0.69 

0.05 

0.33 

0.01 

0.25 

0.04 

0.29 

0.06 

3E-18 

0.64 

0.02 

1.16 

0.58 

-0.67 

0.10 

0.25 

0.04 

0.27 

0.09 

0.37 

0.0 

L08 

rest-UV 













30 

0.61 

0.04 

2.81 

0.51 

-0.84 

0.06 

0.10 

0.04 

0.39 

0.09 

0.37 

0.04 

29 

0.60 

0.03 

2.58 

0.55 

-1.02 

0.12 

0.05 

0.04 

0.68 

0.12 

0.25 

0.11 

28 

0.55 

0.03 

2.12 

0.61 

-1.07 

0.22 

0.01 

0.02 

0.72 

0.13 

0.54 

0.10 

rest-optical 

29 

0.59 

0.03 

2.71 

0.03 

-1.21 

0.06 

0.37 

0.00 

0.07 

0.03 

0.20 

0.01 

28 

0.59 

0.03 

2.69 

0.07 

-1.24 

0.09 

0.16 

0.02 

0.13 

0.03 

0.16 

0.02 

27 

0.60 

0.02 

2.63 

0.11 

-1.25 

0.07 

0.08 

0.01 

0.25 

0.04 

0.16 

0.05 

26 

0.59 

0.02 

2.50 

0.20 

-1.26 

0.09 

0.04 

0.01 

0.47 

0.07 

0.17 

0.07 

25 

0.54 

0.02 

2.50 

0.20 

-1.17 

0.09 

0.04 

0.01 

0.62 

0.10 

0.24 

0.05 

hya 

5E-19 

0.55 

0.00 

2.77 

0.35 

-1.75 

0.12 

0.29 

0.03 

0.62 

0.03 

0.15 

0.07 

3E-18 

0.45 

0.01 

2.06 

0.43 

-1.33 

0.25 

0.10 

0.04 

0.74 

0.11 

0.37 

0.09 

8E-18 

0.41 

0.11 

1.81 

2.61 

-0.75 

0.81 

0.04 

0.27 

0.74 

0.82 

0.54 

0.70 

Ua 2E-19 

0.60 

0.01 

1.67 

0.04 

-0.91 

0.04 

0.36 

0.01 

0.23 

0.02 

0.16 

0.06 

6E-19 

0.57 

0.01 

1.57 

0.12 

-0.87 

0.02 

0.33 

0.01 

0.24 

0.01 

0.13 

0.02 

lE-18 

0.56 

0.01 

1.63 

0.17 

-0.88 

0.03 

0.28 

0.01 

0.29 

0.03 

0.11 

0.02 

3E-18 

0.50 

0.01 

1.67 

0.22 

-1.00 

0.13 

0.11 

0.02 

0.60 

0.13 

0.47 

0.11 


Notes, continuation of Tab. [9] 
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L. Guaita: LARS at low and high z 

Table 11: Morphological parameters of high-z simulated LARS galaxies 


depth 

QrPlO 


QminA 

O-Cn^nA 

M20 

O-M20 

A 

crA 

S 

a-S 

ell 

crell 

L09 

rest-UV 













30 

0.70 

0.02 

3.59 

0.63 

-1.23 

0.10 

0.18 

0.05 

0.46 

0.06 

0.79 

0.08 

29 

0.61 

0.03 

1.54 

1.09 

-0.90 

0.22 

0.19 

0.08 

0.46 

0.10 

0.57 

0.17 

28 

0.69 

0.05 

2.46 

0.56 

-1.74 

0.65 

0.04 

0.03 

0.51 

0.25 

0.38 

0.17 

rest-optical 

29 

0.67 

0.01 

2.85 

0.05 

-0.98 

0.01 

0.36 

0.00 

0.33 

0.00 

0.81 

0.00 

28 

0.65 

0.01 

2.84 

0.12 

-0.98 

0.01 

0.21 

0.01 

0.36 

0.01 

0.83 

0.00 

27 

0.62 

0.06 

2.23 

0.71 

-1.01 

0.26 

0.15 

0.05 

0.50 

0.13 

0.57 

0.22 

Lyor 

5E-19 

0.47 

0.00 

0.82 

0.10 

-0.84 

0.01 

0.28 

0.07 

0.57 

0.02 

0.48 

0.02 

Ua 

2E-19 

0.70 

0.02 

4.38 

0.27 

-1.46 

0.14 

0.23 

0.03 

0.33 

0.12 

0.67 

0.03 

6E-19 

0.67 

0.01 

4.26 

0.29 

-1.49 

0.23 

0.19 

0.01 

0.40 

0.04 

0.79 

0.00 

lE-18 

0.63 

0.03 

4.13 

0.29 

-1.49 

0.50 

0.16 

0.02 

0.46 

0.09 

0.71 

0.15 

3E-18 

0.58 

0.03 

3.69 

0.61 

-2.04 

0.84 

0.11 

0.03 

0.49 

0.16 

0.59 

0.12 

lE-17 

0.56 

0.04 

2.96 

0.60 

-1.81 

0.44 

0.06 

0.04 

0.58 

0.20 

0.36 

0.14 

LIO 

rest-UV 













30 

0.59 

0.02 

2.11 

0.25 

-1.27 

0.14 

0.14 

0.02 

0.21 

0.04 

0.54 

0.02 

29 

0.57 

0.03 

2.33 

0.41 

-1.30 

0.13 

0.07 

0.02 

0.37 

0.09 

0.60 

0.05 

28 

0.53 

0.04 

2.27 

0.60 

-1.31 

0.23 

0.02 

0.02 

0.57 

0.15 

0.62 

0.09 

rest-optical 

29 

0.57 

0.01 

3.34 

0.26 

-1.97 

0.04 

0.12 

0.00 

0.09 

0.00 

0.55 

0.02 

28 

0.59 

0.01 

3.35 

0.16 

-2.04 

0.04 

0.08 

0.01 

0.20 

0.02 

0.50 

0.02 

27 

0.59 

0.01 

3.18 

0.21 

-1.98 

0.09 

0.04 

0.01 

0.41 

0.05 

0.52 

0.03 

26 

0.57 

0.01 

2.95 

0.34 

-1.92 

0.10 

0.02 

0.01 

0.62 

0.08 

0.53 

0.07 

hya 

5E-19 

0.45 

0.06 

1.08 

0.40 

-0.82 

0.66 

0.17 

0.13 

0.74 

0.18 

0.41 

0.11 

Ua 

2E-19 

0.58 

0.01 

2.38 

0.70 

-1.06 

0.14 

0.24 

0.02 

0.30 

0.06 

0.59 

0.03 

6E-19 

0.58 

0.01 

2.38 

0.70 

-1.06 

0.14 

0.24 

0.02 

0.30 

0.06 

0.59 

0.03 

lE-18 

0.55 

0.02 

2.56 

0.68 

-1.05 

0.15 

0.18 

0.03 

0.35 

0.09 

0.65 

0.03 

Lll 

rest-UV 













30 

0.57 

0.02 

3.20 

0.10 

-1.21 

0.08 

0.30 

0.04 

0.39 

0.03 

0.88 

0.08 

29 

0.56 

0.04 

2.58 

0.53 

-1.11 

0.18 

0.06 

0.05 

0.47 

0.07 

0.88 

0.09 

28 

0.57 

0.04 

2.53 

0.74 

-0.70 

0.18 

0.04 

0.03 

0.69 

0.10 

0.62 

0.15 

rest-optical 

29 

0.51 

0.01 

2.95 

0.46 

-1.23 

0.10 

0.28 

0.07 

0.26 

0.13 

0.82 

0.13 

28 

0.51 

0.02 

3.03 

0.08 

-1.33 

0.03 

0.18 

0.06 

0.28 

0.01 

0.87 

0.00 

27 

0.49 

0.02 

2.94 

0.17 

-1.31 

0.07 

0.05 

0.05 

0.43 

0.05 

0.87 

0.05 

26 

0.51 

0.02 

2.60 

0.37 

-0.89 

0.12 

0.01 

0.01 

0.64 

0.07 

0.75 

0.12 

hya 

5E-19 

0.44 

0.34 

1.62 

1.54 

-0.73 

1.69 

0.43 

0.32 

0.72 

0.87 

0.63 

0.98 

3E-18 

0.39 

0.05 

1.81 

0.54 

-0.65 

0.16 

0.29 

0.07 

0.63 

0.25 

0.64 

0.07 

Ha 

2E-19 

0.58 

0.02 

2.61 

0.08 

-0.89 

0.10 

0.23 

0.01 

0.37 

0.03 

0.88 

0.01 

6E-19 

0.52 

0.01 

2.39 

0.16 

-1.11 

0.14 

0.13 

0.02 

0.37 

0.11 

0.81 

0.18 

lE-18 

0.48 

0.10 

2.58 

0.46 

-0.67 

0.36 

0.10 

0.12 

0.73 

0.29 

0.44 

0.64 

3E-18 

0.44 

0.01 

2.28 

0.40 

-1.01 

0.11 

0.03 

0.02 

0.63 

0.10 

0.73 

0.13 


Notes, continuation of Tab. 
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L. Guaita: LARS at low and high z 

Table 12: Morphological parameters of high-z simulated LARS galaxies 


depth 

QrPlO 


QminA 

O-Cn^nA 

M20 

O-M20 

A 

crA 

S 

a-S 

ell 

crell 

L12 

rest-UV 

30 

0.69 

0.00 

3.33 

0.23 

-1.43 

0.03 

0.37 

0.01 

0.10 

0.02 

0.38 

0.01 

29 

0.69 

0.01 

3.32 

0.24 

-1.44 

0.04 

0.25 

0.03 

0.11 

0.01 

0.41 

0.01 

28 

0.70 

0.01 

3.29 

0.30 

-1.47 

0.22 

0.16 

0.04 

0.18 

0.03 

0.45 

0.02 

27 

0.69 

0.03 

3.11 

0.53 

-1.84 

0.41 

0.11 

0.04 

0.39 

0.19 

0.27 

0.17 

rest-optical 

29 

0.65 

0.01 

3.09 

0.22 

-1.71 

0.15 

0.37 

0.07 

0.15 

0.01 

0.49 

0.01 

28 

0.65 

0.01 

3.20 

0.30 

-1.65 

0.10 

0.25 

0.01 

0.15 

0.01 

0.52 

0.01 

27 

0.66 

0.01 

3.35 

0.39 

-1.80 

0.18 

0.14 

0.01 

0.20 

0.01 

0.54 

0.01 

26 

0.68 

0.05 

3.21 

0.60 

-1.54 

0.48 

0.07 

0.04 

0.41 

0.12 

0.54 

0.09 

Lyor 

5E-19 

0.62 

0.01 

1.27 

0.11 

-2.26 

0.25 

0.37 

0.03 

0.44 

0.10 

0.27 

0.03 

3E-18 

0.54 

0.06 

1.38 

0.30 

-1.81 

0.87 

0.13 

0.07 

0.75 

0.15 

0.47 

0.13 

Ua 

2E-19 

0.61 

0.05 

2.99 

0.11 

-1.54 

0.22 

0.34 

0.05 

0.19 

0.14 

0.39 

0.06 

6E-19 

0.62 

0.01 

3.01 

0.00 

-1.66 

0.16 

0.28 

0.00 

0.19 

0.06 

0.40 

0.02 

IE-18 

0.62 

0.01 

3.01 

0.00 

-1.83 

0.20 

0.26 

0.01 

0.17 

0.05 

0.41 

0.02 

3E-18 

0.60 

0.01 

2.93 

0.29 

-1.60 

0.30 

0.17 

0.01 

0.23 

0.03 

0.44 

0.03 

L13 

rest-UV 

30 

0.72 

0.03 

1.60 

0.36 

-1.24 

0.34 

0.27 

0.01 

0.17 

0.12 

0.37 

0.11 

29 

0.71 

0.03 

2.30 

0.73 

-1.35 

0.24 

0.23 

0.03 

0.23 

0.13 

0.38 

0.08 

28 

0.65 

0.03 

2.36 

0.65 

-1.37 

0.19 

0.18 

0.03 

0.20 

0.12 

0.15 

0.08 

rest-optical 

29 

0.64 

0.02 

3.64 

0.03 

-1.32 

0.22 

0.31 

0.00 

0.24 

0.15 

0.47 

0.06 

28 

0.62 

0.01 

3.52 

0.08 

-1.45 

0.06 

0.21 

0.01 

0.23 

0.01 

0.53 

0.02 

27 

0.49 

0.05 

3.31 

0.11 

-0.98 

0.16 

0.16 

0.01 

0.49 

0.15 

0.53 

0.11 

Lya 

5E-19 

0.53 

0.03 

2.10 

1.28 

-0.52 

0.12 

0.40 

0.07 

0.61 

0.09 

0.40 

0.12 

3E-18 

0.48 

0.02 

1.74 

0.50 

-1.18 

0.32 

0.28 

0.07 

0.55 

0.09 

0.24 

0.10 

8E-18 

0.44 

0.02 

1.63 

0.56 

-1.22 

0.25 

0.12 

0.06 

0.72 

0.13 

0.40 

0.11 

Ua 

2E-19 

0.67 

0.01 

1.99 

0.00 

-1.90 

0.09 

0.16 

0.00 

0.17 

0.03 

0.23 

0.07 

6E-19 

0.66 

0.01 

1.99 

0.00 

-1.87 

0.10 

0.16 

0.01 

0.24 

0.03 

0.25 

0.07 

lE-18 

0.65 

0.02 

1.99 

0.00 

-1.93 

0.43 

0.14 

0.01 

0.30 

0.11 

0.28 

0.17 

3E-18 

0.60 

0.02 

1.80 

0.30 

-1.91 

0.32 

0.09 

0.02 

0.44 

0.15 

0.34 

0.16 

L14 

rest-UV 

30 

0.64 

0.00 

2.39 

0.00 

-0.99 

0.04 

0.36 

0.00 

0.04 

0.02 

0.41 

0.02 

29 

0.64 

0.01 

2.39 

0.00 

-0.77 

0.11 

0.32 

0.01 

0.11 

0.07 

0.52 

0.03 

28 

0.65 

0.05 

2.38 

0.09 

-0.94 

0.20 

0.29 

0.05 

0.08 

0.10 

0.47 

0.15 

rest-optical 

29 

0.78 

0.00 

2.39 

0.00 

-1.06 

0.09 

0.39 

0.00 

0.19 

0.01 

0.09 

0.02 

28 

0.78 

0.01 

2.39 

0.00 

-1.17 

0.12 

0.29 

0.02 

0.27 

0.04 

0.22 

0.05 

27 

0.67 

0.13 

1.97 

0.51 

-0.93 

0.34 

0.28 

0.11 

0.36 

0.18 

0.50 

0.12 

Lya 

5E-19 

0.61 

0.03 

2.13 

0.00 

-1.47 

0.07 

0.29 

0.00 

0.16 

0.01 

0.17 

0.01 

3E-18 

0.61 

0.02 

2.13 

0.01 

-1.51 

0.09 

0.24 

0.01 

0.20 

0.03 

0.24 

0.04 

8E-18 

0.56 

0.03 

2.12 

0.04 

-1.48 

0.09 

0.19 

0.01 

0.28 

0.07 

0.21 

0.04 

2E-17 

0.50 

0.02 

2.04 

0.16 

-1.44 

0.09 

0.12 

0.02 

0.40 

0.11 

0.15 

0.05 

5E-17 

0.43 

0.02 

1.99 

0.23 

-1.33 

0.12 

0.05 

0.03 

0.60 

0.11 

0.23 

0.09 

Ha 

2E-19 

0.55 

0.01 

1.50 

0.00 

-1.07 

0.03 

0.20 

0.00 

0.03 

0.01 

0.39 

0.01 

6E-19 

0.57 

0.01 

1.50 

0.00 

-1.08 

0.01 

0.21 

0.00 

0.03 

0.00 

0.42 

0.01 

lE-18 

0.58 

0.01 

1.50 

0.00 

-1.09 

0.01 

0.18 

0.00 

0.04 

0.00 

0.45 

0.02 

3E-18 

0.65 

0.07 

1.63 

0.21 

-0.91 

0.30 

0.19 

0.09 

0.23 

0.18 

0.43 

0.13 

lE-17 

0.52 

0.03 

1.48 

0.12 

-1.35 

0.12 

0.10 

0.02 

0.14 

0.07 

0.19 

0.08 


Notes, continuation of Tab. [9] 
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L. Guaita: LARS at low and high z 


Appendix A: Non-parametric measurements 

In this section we explain the way we measured galaxy sizes 
and non-parametric quantities in detail, and show the results for 
galaxies with known-profiles. 

In Fig. |A.l| we summarize the equations adopted in this anal¬ 
ysis and first introduced by|Conselice et al.|(|2000|) and|Lotz et al.| 
( |2004l ). 

The code we developed makes a basic use of the 
ELLIPSE task in iraf.stsdas.isophote and the PHOT task 
in iraf.digiphot.apph ot. We first ran the So urce Extractor 
(SExtractor) software ( [Bertin & Arnouts||1996| ) on one galaxy 
image. This provided the centroid and the elliptical aperture 
containing the entire galaxy. We adopted configuration param¬ 
eters like in [Bond et al.| ( |200^ (DETECT.THRESH =1.65 
and DEBLENDJV[INCONT=l). We followed the choice of 
DETECT_MINAREA=30 for the high-z simulations. Those pa¬ 
rameters were optimized to provide significant morphological 
measurements in deep HST-band observations. A larger value of 
contiguous pixels was adopted to prevent SExtractor from break¬ 
ing up the clumpy, resolved, original z ~ 0 LARS galaxies into 
smaller fragments. 

We adopted SExtractor centroid, orientation angle, and el- 
lipticity as the fixed ELLIPSE parameters and the SEx AUTO 
photometry semi-major axis as the reference semi-major axis 
length (smaO). We then measured flux within ellipses by vary¬ 
ing the semi-major axis (sma). The task was able to fit elliptical 
isophotes at a pre-defined, fixed sma, and works better for well- 
defined galaxy profiles. As LARS galaxies are irregular, a better 
convergence of the task was performed by fixing the ellipse ori¬ 
entation and ellipticity. 

The ELLIPSE task outputs surface brightness (I(r)) and inte¬ 
grated flux (E) within every sma (r^, r^+i, ...). W e used the given 
surfac e brightness to derive the Petrosian ratio ( [Bershady et al.| 
|2000[ Tj =I(r)/<I(r)>) as a function of sma and the integrated 
flux to estimate r 2 o, Uo. rgo, the radii containing 20%, 50%, 
and 80% of the total source flux. The Petrosian semi-major axis 
(rP20^^^) corresponds to the sma where p = 0.2. We defined an 
elliptical concentration proportional to rgo/Uo- Applying a 
smoothing kernel with width equal to rP20/5, we also estimated 
a smoothed-image Petrosian radius (rP20S, Lotz et al.|2004| ). 

The PHOT task outputs fiuxes integrated within circular 
apertures. We derived the corresponding I(r) and estimated 
the circular Petrosian radius (rP20‘^*^‘^), the circular r^Q^, and 

) was then proportional to 


rgQ^. The circular concentration (C 

^circ fyCirc 
^80 /^20 • 

We also defined the signal-to-noise (SN) per pixel as 


SN 


pixel 




and calculated the Petrosian radius (rP20'^^"^), and 
corresponding to the minimum of asymmetry. The point of the 
galaxy of minimum asymmetry is generally close to its brightest 
pixel, but not necessarily to the SExtractor centroid. The con¬ 
centration at mi nimum asymmetry was then calculated 


from the ry^/r^^^^ ratio (Conselice et a l.|2000 
Jiang et al.|2013[|Holwerda et al.|2(JI4p 


Lotz et al. 


2006 


Bershady et al.| ( |2000| ) defined C by measuring photometry 


inside circular apertures. They estimated that C could be under¬ 
estimated up to 30%, in the case of an ellipticity of ~ 0.75 and 
circular aperture, but that it was within 10-15% for early-type 
galaxies. was the quantity we mainly used throughout the 
paper. 

A galaxy was assumed to be composed of a certain number 
of pixels, constituting a segmentation map. The non-parametric 
measurements and SN estimation were performed counting the 
fiux of pixels belonging to that segmentation map. We defined 
the segmentation map in two ways, one (the fixed-size segmen¬ 
tation map) i s an ellipse with sma = rP20 and orient ation given 
by SE x (see Scarlata et al.||200^ ; one described in Lotz et al. 
( |2004| ), where the pixels belonging to the segmentation map have 
surface brightness larger than the value at rP20S. The fixed-size 
map was mainly concentrated in the central part of the galaxy, 
the second one could contain bright pixels in the galaxy out¬ 
skirts. 

The Gini coefficient (G) a nd M 2Q were also calculated by 
following the equations in Eig. |A.1[ X, (i = 1 to n) correspond 
to the pixel values, sorted in increasing order, and X the aver¬ 
age pixel value, within the chosen segmentation map. f/ (i = 1 
to n) correspond to the pixel values within the chosen segmenta¬ 
tion map, but sorted in decreasing order; x^ and y^ are the pixels 
corresponding to SEx centroid. A first value of G (G^^^^) was 
estimated within the fir st segmentation map dLotz et al.||2()04| 


2006[ Jiang et al.|2013t |Holwerda et al.|2014| i. A second value 

(G^^^ ) was estimated within the second segmentation map. 

The latter is sensitive to multiple knots in a galaxy full of struc¬ 
tures. Within the fixed-sized segmentation map we measured 


where CTsky is the standard deviation of means measured in 
more than three boxes around the galaxy, S/ the signal, and n 
the number of pixels belonging to a galaxy. The total SN of the 
galaxy was then obtained by multiplying the SNpixei by ^/n. If 
> | 5 ,|, 

The asymmetry (A) was calculated as the minimum value of 
the normalized difference between the galaxy image (Iq) and the 
same one rotated by 1 80^ (Iiso)- We adopted t he background (B) 
correction advised by Conselice et al.| ( |2000[ ) for low SN galax- 
ies. Based on this definition, 0 < A < 1 (e.g., [Scarlata et al.| 
2007t|Aguirre et al.|2013[|Law et al.|2012| ). After calculating the 


position of minimum asymmetry, we ran PHOT on that position 


M2Q ([Scarlata et al.|2007 pTang et al.|2013t|Aguirre et al.|2013[ 

[Holwerda et al.|2Q14[) and S. 

The dumpiness (S) was defined by [Conselice[ ( [2003' ) as the 
normalized difference between the galaxy image (I) and the 
smoothed image (Io. 3 xrP 205 where the smoothing kernel sigma 
was 0.3X rP20). The pixels belonging to the galaxy image are 
the ones within 0.3 and 1.5 times the rP20, i.e. we excluded the 
very central pixels, which are often unresolved. 

We tested our code on the frames showed in Pig. [A.2| The 
code calculations are shown in Pig. [A.3[ Pirst of all we analyti¬ 
cally calculated the theoretical (THEO in the figure) Petrosian 
ratio and radius of an exponential and deVacouleur profiles. 
Then, we ran SExtractor and calculated the non-parametric mea¬ 
surements described above. 

The results are listed in Tables IajI and [A.2[ The rP20^^^ size 
better recovers the analytical value of an exponential profile. It 
also well recovers the value in the case of added noise and of a 
deVacouleur profile. Por this reason, we adopted rP20^^^ as the 
main size estimator throughout the paper. By comparing with 
the an alyti cal solution and th e estimations by [Bershady et al. 
( [200Q[ ) and [Lotz et al.[ ( |2006| ), we noticed that we could under¬ 
estimate of a deVacouleur profile up to 30%, of an 
exponential(deVacouleur) profile up to 10(20)%, and overesti¬ 
mate up to 10(5)% for an exponential(deVacouleur) 

profile. M20 is well-recovered for all the profiles within 3%. 

Also, we tested our code on the public galaxy stamps from 
COSMOS and compared the results with the ones obtained with 
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n(r)= rP20 1 n(rP20)=0.2 Consdiceiooo.ioo} 

(/(r)) Ar,^,-Ar, Fr, 

i J i 21^0 J 


C = 5 X log(r80%/''2(y>/.) 


G = - ' f (2i n \)Xi 

Xj within the segnnentation nnap 

M 20 = while ^75 <0.27;ot 

fj within the segnnentation map 

[(XI - x,f + (y,- - y,f] 


1.5xrP20/T J \ 

g _ ^ 0 ^5=0.3xrP20/ 

0.3X/-P20 ^ 



Fig. A.l: Basic equations adopted to calculate galaxy size and non- 
parametric measurements of the LARS original and high-z simulated 
galaxies. 





Fig. A.2: Simulated profiles used as a test of our code performance. 
From the upper left to the lower right panels: symmetric exponential 
profile, exponential profile with ellipticity equal to 0.5 and position an¬ 
gle 0^, exponential profile with ellipticity equal to 0.5 and position angle 
45^, deVacouleur profile, asymmetric exponential profile after adding a 
noise corre spo nding to a lOcr detection limit 7 times fainter (noise in 
Tables [AT] and |A.2] ) and 3 times fainter than the fake galaxy fiux. They 
were generated by running the MKOBJECT task in iraf.artdata and 
noise was added by running the MKNOISE task in iraf.artdata. 



the ZEST software ( [Scarlata et al.||2007] ). When assuming the 
brightest pixel as the centre of a galaxy, we recovered val¬ 
ues in more than 80% of the cases. As the fixed-size s egmenta¬ 
tion map is a better choice for redshift comparisons ( Scarlata] 
et al. 2007), we tended to prefer rather than 


throughout the paper. 


owing pixel resampling and added noise, we defined A(param), 


A{param) = 


param(z ~ 2, noise) - param(original) 
param(original) 


(A.1) 


which represented the difference between the measurement 
of a parameter in the simulated and in the original image, 
normalized to the value in the original image. By definition, 
A(param) tended to be zero when the two measurements were 
very similar, tended to be equal to -1 when the measurement in 
the simulated image was much smaller than in the original one, 
and equal to 1 when the measurement in the simulated image 
was twice that in the original one. We calculated A (param) for 
every gala xy an d estimated the mean and the standard deviation. 

In Fig. A.4 we show the results. A, and M20 

were all preserved after pixel resampling in the continuum im¬ 
ages. Owing survey depth, A decreased. 

[Huertas-Company et al.| ( |2Q14] ) studied the variation of A, G, 
and M20 due to resampling the rest-frame optical of local galax¬ 
ies to z > 1. They found that A can increase up to 50%, G up 
to 10%, while M20 can decrease up to 10%. These trends were 
observed to be more pronounced for early-type galaxies; Lotz 
et al. I poo's] ) resampled the rest-frame UV of a few sources from 
z = 1.5 to z = 4, findin g that G and M20 were preserved, con¬ 
sistent with our results. Qverzier et al. ( ]20101 ) investigated the 
change of A, C, G, and M20 when resampling z = 0.2 Lyman 
break analogues to z = 2. They gave an R-value scale, where 
|R| ~ 0 when the difference between the median parameter mea¬ 
sured in the resampled and in the original image was small and 
|R| ~ 1 when was comparable to the sample scatter. They found 
= 0.02, = -0.34, and R^^^ = -0.43 in the 

rest-frame UV and = -0.46, R^ = -0.52, R^ = -1.1, and 


R^^^ = -0.26 in the rest-frame optical. Recently, Petty et al. 
( 2014] ) explored the variations of G and M20 of local LIRGS 
when simulated to be at z = 0.5, 1.5, 2, 3 and in a survey with 
the HUDF depth. Some of the galaxies in their sample, charac¬ 
terized by clumps and filamentary structures (typical of merg¬ 
ing systems) in high-resolution HST images, tended to appear 
as disk-like galaxies by z = 2. Some others maintained merging 
systems’ morphology. All these trends are in agreement with our 
findings, but they also tell us that the variations do not follow a 
specific trend for irregular galaxies. 

In the Lya images, A tended to decrease and could in¬ 
crease due to pixel resampling. Owing survey depth, all the four 
parameters tended to decrease and A significantly decreased. 

Another way to study the change of morphological param¬ 
eters is to loo k at p airs of them, like G^^^^ vs M20 and A vs 


-^minA 


In Fig. 


A.5 


we show the diagrams of two-parameters. 


measured in the rest-frame optical, UV, and Lya images. Our 
analysis showed that the G and M20 estimations were not af¬ 
fected by resolution and survey depth. Therefore, they were use¬ 
ful for comparisons at different redshifts and different survey 
depths throughout the paper. 


A.1. Comparison between original and simulated galaxy 
morphological parameters 

Following the method described in the previous section, we esti¬ 
mated sizes and morphological parameters in the continuum and 
line images. We adopted the same method in the case a galaxy 
was simulated to be at z ~ 2 and observed in different depth 
surveys. To quantify the morphological parameter variations 


34 






































L. Guaita: LARS at low and high z 



Fig. A.3: From the left to the right panels surface brightness, Petrosian ratio, and integrated flux. The upperQower) row is obtained by using 
elliptical(circular) aperture and the iraf ELLIPSE(PHOT) task. Eor an effective radius. Re = 10 pixels, we estimated photometry in the case of a 
symmetric exponential profile (solid line), exponential profile with ellipticity equal to 0.5 and either position angle 0 (dashed line) or 45 (dotted 
line) degrees, and of a deVacouleur profile (triangles). The case of adding noise is represented by filled circles with and without errorbars. The 
analytical Petrosian ratios for exponential (THEOexp) and deVacouleur (THEOdeVauc) profiles are presented in the second-column panels as blue 
and cyan solid curves respectively. 


Table A.l: Sizes from known profile galaxies of Re=10 pixels 


profile 

rP20^'' 

rP20"'"" 

rP20^*^ 

r20 

rso 

rso 

kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

kpc 

exp 

ellO 

22.0 

21.0 

19.0 

5.5 

10.0 

17.5 

ell0.5,pa0 

21.5 

16.0 

15.0 

5.0 

10.0 

17.5 

ell0.5,pa45 

21.5 

16.0 

15.0 

5.0 

10.0 

17.0 

noise 

19.0:^ 

16.0^^ 

23.0 

5.0 

10.0 

16.0 

deVauc 

20.5 

17.0 

11.0 

2.5 

7.5 

17.0 


Notes. The analytical values for the ratio of rP20 and effective radius is 2.16 for an 

ITO - 


exponential and 1.82 for a deVacouleur profile ( [Bershady et al.| 


Table A.2: Morphology measurements from known profile galaxies of Re = 10 pixels 


profile 

QrPlO 

QSB-rP2i)S 

M20 

ell 

Qcirc 

QeU 

O^minA 

A 

S 

exp 










ellO 

0.4888 

0.5229 

-1.752 

0.05 

2.657 

2.513 

2.386 

0.0 

0.016 

ell0.5,pa0 

0.4923 

0.5326 

-1.836 

0.5 

2.559 

2.720 

2.559 

0.0 

0.036 

ell0.5,pa45 

0.4874 

0.5333 

-1.796 

0.5 

2.559 

2.657 

2.559 

0.0 

0.034 

noise 

0.4914 

0.5164 

-1.164 

0.114 

2.236 

2.526 

2.236 

-0.0135 

0.393 

deVauc 











0.6183 

0.6271 

-2.398 

0.05 

3.495 

4.162 

3.495 

0.0 

0.025 


Notes . The analytical values of is 2.7 and 5.2 for an exponential and a deVacouleur profile respectively ( |Ber shady et al. 2000| ). Lotz 
( [2006| ) estimated =0.473, M20=-1.8, C™”^=2.71 for a noise-free image of an exponential profile and =0.6, M20=-2.47, 


et al. 


=4.34 for a noise-free image of a deVacouleur profile (half-light-radius equal to 600 resolution elements). 
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rest-frame optical survey depth rest-frame optical survey depth 



rest-frame UV survey depth rest-frame UV survey depth 


Lvo line 





10 ^ 10 " 
Fdya) survey depth x lE-18 


Fig. A.4: Mean A(param) as defined in eq. A.l vs simulated survey depth, param is A {upper left), {upper right), {lower left), and 


M20 {lower right) for the rest-frame optical, UV, and Lycr simulated images. The deepest simulated survey is represented by the most right data 
point in each panel. The error bars represent the standard deviation among all the high-z simulated LARS galaxies. In the case fewer-than-twelve 
galaxies are detected at a specific depth, the error bars is increased proportionally to the number of undetected galaxies. The numbers on the top 
right corner of each panel represent the Spearman test coefficient, r, and probability, p, of uncorrelated datasets, assuming that the depth vector 
indicates the deepest survey on the right of the x-axis. 
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L. Guaita: LARS at low and high z 



Fig. A.5: vs M20 (first column) and A vs (second column) measured in the rest-frame optical {upper), UV {middle), and Lycr {lower) 

images. Small and big diamonds with error bars indicate the parameters of the original and (deepest survey) high-z simulated LARS galaxies. The 
dashed and dotted-dashed lines indicate the separation between merging system and normal galaxies’ parameter space. Owing pixel resampling 
the galaxies stay in the same region of the vs M20 diagram. 
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Appendix B: LARS galaxies simulated at high 
redshift in a deep and shallow survey 


frame UV and optical continua, Lya and Ha lines. The profiles 
are normalized to 2 kpc to compare continuum and line profiles. 


In this appendix, we show the continuum and line maps of 


present rest-frame UV, Lya, rest-frame optical, and Ha images 
in the deepest survey depth probed here, together with SEx de¬ 
tection apertures. In Figs. |B.4[|B.5| and |B.6| we show the same 
stamps but for a shallower simulated survey. 


the high-z simulated LARS galaxies. Figures [BT]|B.2[ and B3 


Appendix C: Surface brightness profiles of original 
and high-z simulated LARS galaxies 


We present here the surface brightness profiles of eleven LARS 
galaxies studied in this paper. The profiles of LOl were shown 
in Fig.[^ Every figure in this appendix shows four panels, rest- 
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Fig. B.l: Rest-frame UV (first column), Lycr (second column), rest-frame optical (third column), and Her (fourth column panels) images for z2L01, 
z2L02, z2L03, z2L05, and z2L07 in the deepest survey, probed here. Every panel is 20x17 kpc wide. The reddish ellipses indicate SEx apertures 
(dashed curves indicate fiagged sources according to SEx convention), corresponding to the assumed detection parameters: DETECT THRESH 
= 1.65, DETECT MINAREA=30, and DEBLEND MINCONT=l from Bond et al. (2009). The colour scaling is logarithmic and chosen to show a 
visually consistent background noise. 
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L. Guaita: LARS at low and high z 



Fig. B.2: As Fig. B.l for z2L08, z2L09, z2L10, z2Lll, and z2L12. In the case of z2L09 and z2Lll the size is 27x27 kpc and 33x30 kpc respectively 
to fit their elongated shapes. For these galaxies SEx apertures in Lya happen to be outside the shown region. 
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F(Lya)=5E-19 
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F(Ha)=6E-19 
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Fig. B.3: As Fig. B.l for z2L13 and z2L14. For z2L13 SEx aperture in Lycr happens to be outside the shown region. Also, two main sources 
connected by undetectable {mrest-uv >30 and F(Hcr) < 6E-19 erg sec"^ cm"^) surface brightness structures are seen in UV continuum and Her. By 
using the chosen detection parameters, SEx found two sources as separated. As in Lycr and optical continuum the photometric measurements are 
done centring the aperture close to the right clump, we locate the aperture on that one for the photometry in UV and Her as well. 
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Fig. B.4: As B .l, but for a shallow simulated survey. The dashed-line apertures indicate that some detected source is blended to another ( [Bertin & 


Arnouts 
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Fig. B.5: As Fig. B.2 but for a shallow simulated survey. z2L09 and z2L10 are not detected in Lycr at a lOcr detection limit of 3E-18 erg sec 
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F(Lya)=3E-18 
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Fig. B.6: As Fig. |B.3| but for a shallow simulated survey. Only the bright right clump of z2L13 is detected in UV at a lOcr detection limit of 
^rest-uv - 29. Its Lyof emission is localized around the lower side of the right clump. In Her the two overlapping sources become well separated. 




Fig. C.l: Normalized surface brightness profiles of L02. Black points with error bars correspond to the surface brightness profile of the original 
LARS images in the rest-frame UV, optical, and Lycr as explained in the text (Fig.[^. The red squares represent the profile affected by background 
noise, for a certain shallow survey. 
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sma [kpc] 





Fig. C.2: Same colour coding as in Fig. CJ_ but for z2L03 


in magenta(red) symbols with error bars, as throughout the 


and z2L05. The deepest and 
paper. 


intermediate depth 


surveys are shown for z2L03(z2L05) 








Fig. C.3: Same colour coding as in Fig. C.l but for z2L07 
in blue(light green) symbols with error bars. 


and z2L08. The deepest and intermediate depth 


surveys are shown for z2L07(z2L08) 









Fig. C.4: Same colour coding as in Fig. |C.l| but for z2L09 and z2L10. The deepest and intermediate depth surveys are shown for z2L09(z2L10) 
in black(dark red) symbols with error bars. 
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L. Guaita: LARS at low and high z 








Fig. C.5: Same colour coding as in Fig. CJ_ but for z2Ll 1 


in dark green(orange) symbols with error bars. 


and z2L12. The deepest and intermediate depth surveys are shown for z2Lll(z2L12) 



■ L13original 

1 L13 m(opticalcont)=29.0 
L13 m=28.0 

1 L13 m=27.0 









Fig. C.6: Same colour 
in grey(pink) symbols 


coding as in Fig. C.l but for z2L13 
with error bars. 


and z2L14. The deepest and intermediate depth surveys are shown for z2L13(z2L14) 
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